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ABSTRACT 


Polyphase  continuous  waveform  (CW)  radar  systems  often  use  the  popular  Frank  code 
and  P4  code  due  to  their  linear  time-frequency  characteristics  as  well  as  their  low 
periodic  ambiguity  sidelobes.  The  phase  relationship  of  the  Frank  code  corresponds  to  a 
sawtooth  folding  waveform.  The  phase  relationship  of  the  P4  code  is  symmetrical  with  a 
parabolic  distribution.  The  radar  system’s  unambiguous  target  detection  range  is  limited 
by  the  number  of  subcodes  within  the  code  period  (code  length).  Increasing  the  code 
length  to  extend  the  unambiguous  range  results  in  a  larger  range-Doppler  correlation 
matrix  processor  in  the  receiver,  a  longer  compression  time  and  an  increase  in  the 
receiver’s  bulk  memory  requirements.  In  addition,  the  entire  code  period  may  not  be 
returned  from  the  target  due  to  a  limited  time-on-target  resulting  in  significant  correlation 
loss.  To  significantly  extend  the  unambiguous  range  beyond  a  single  code  period,  this 
thesis  explores  the  relationship  between  the  polyphase  codes  (Frank  and  P4) 
and  the  number  theoretic  transforms  (NTT)  where  the  residues  exhibit  the  same 
distribution  as  the  polyphase  values.  The  unambiguous  range  is  extended  from  the 
number  of  subcodes  within  a  single  code  period  to  the  dynamic  range  of  the  transform 
without  requiring  a  large  increase  in  correlation  processing.  The  dynamic  range  of  a  NTT 
is  defined  as  the  greatest  length  of  combined  phase  sequences  that  contain  no  ambiguities 
or  repeated  paired  terms.  By  transmitting  N  >  2  coprime  code  periods,  the  unambiguous 
range  can  be  extended  by  considering  the  paired  values  from  each  sequence.  A 
new  Frank  phase  code  formulation  is  derived  as  a  function  of  the  residue  number  system 
(RNS)  where  each  residue  corresponds  to  a  phase  value  within  the  code  period  (modulus) 
sequence.  Based  on  the  symmetrical  distribution  of  the  P4  code,  a  new  phase  code 
expression  is  derived  using  both  the  symmetrical  number  system  (SNS)  and  the  robust 
symmetrical  number  sy stem  (RSNS).  Here  each  phase  value  within  the  code  period 
corresponds  to  a  symmetrical  residue.  MATLAB  simulations  are  used  to  verify  the  new 
expressions  for  the  RNS,  SNS  and  RSNS  phase  codes.  Implementation  considerations  of 
the  new  approach  are  also  addressed. 
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EXECUTIVE  SUMMARY 


The  Frank  phase  code  and  the  P4  phase  code  are  often  used  in  continuous  wave  (CW) 
polyphase  modulation  radar  systems  to  detect  the  target’s  range  due  to  their  linear  time- 
frequency  characteristics  as  well  as  their  low  periodic  ambiguity  sidelobes.  Both  phase 
codes  limit  the  unambiguous  range  Ru  to  within  a  single  code  period  T  =  Nctb ,  where  Nc 

is  the  number  of  subcodes  and  tb  is  the  subcode  period.  To  extend  the  unambiguous 

detection  range,  the  radar  might  simply  increase  the  code  period  of  the  polyphase 
modulation  signal.  However,  the  longer  code  period  requires  a  larger  range-Doppler 
correlation  matrix  processor,  an  increase  in  the  compression  time,  as  well  as  an  increase 
in  the  bulk  memory  in  the  receiver.  Moreover,  a  significant  amount  of  loss  is  incurred  if 
the  total  number  of  code  periods  returned  from  the  target  is  less  than  the  number  of 
reference  code  periods  used  in  the  compression  process  due,  for  example,  to  a  limited 
time-on-target.  A  new  set  of  waveform  architectures  that  can  solve  these  limitations  is 
proposed  in  this  thesis. 

Number  theoretic  transforms  (NTT)  based  on  residue  number  systems  are  used  in 
many  digital  signal  processing  applications  to  increase  the  amount  of  information 
available  from  various  folding  waveforms.  With  a  set  of  N  >  2  moduli,  the  paired  integer 
residue  terms  from  the  transform  sequences  must  be  distinct  (no  ambiguities)  throughout 
the  entire  dynamic  range.  The  dynamic  range  of  a  NTT  is  defined  as  the  greatest  length 
of  combined  sequences  of  integers  that  contain  no  ambiguities  or  repeated  paired  terms. 
Consequently,  the  NTT  residues  within  the  dynamic  range  can  be  used  to  construct  a 
low-power  CW  polyphase  waveform  that  is  useful  for  detecting  targets  beyond  the  single 
code  period  unambiguous  range. 

To  significantly  extend  the  unambiguous  range  beyond  a  single  code  period,  this 
thesis  explores  the  relationship  between  the  polyphase  codes  (Frank  and  P4)  and  the 
NTT  where  the  residues  exhibit  the  same  distribution  as  the  polyphase  values.  The 
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unambiguous  range  is  extended  from  the  number  of  subcodes  within  a  single  code  period 
to  the  dynamic  range  of  the  transform  without  requiring  a  large  increase  in  correlation 
processing. 

A  new  Frank  phase  code  formulation  is  derived  as  a  function  of  the  residue 
number  system  (RNS)  where  each  residue  corresponds  to  a  phase  value  within  the  code 
period  (modulus)  sequence.  Based  on  the  symmetrical  distribution  of  the  P4  code,  a 
new  phase  code  expression  is  derived  using  both  the  symmetrical  number  system  (SNS) 
and  the  robust  symmetrical  number  system  (RSNS).  Here  each  phase  value  within  the 
code  period  corresponds  to  a  symmetrical  residue.  MATLAB  simulations  are  used  to 
verify  the  new  expressions  for  the  RNS,  SNS  and  RSNS  phase  codes.  Implementation 
considerations  of  the  new  approach  are  also  addressed. 

The  relationship  between  the  Frank  phase  code  and  the  RNS  is  developed  based 
on  their  similar  saw-tooth  waveform  distribution.  A  new  equation  using  the  RNS  to 
derive  the  Frank  phase  code  is  formulated.  This  new  phase  code  sequence  is  called  the 
residue-Frank  phase  code.  The  periodic  autocorrelation  function  (PACF)  and  the 
periodic  ambiguity  function  (PAF)  show  that  the  new  residue-Frank  phase  code  formula 
gives  exactly  the  same  results  as  using  the  formula  in  the  Frank  phase  code.  Next,  a  set  of 
N  =  3  RNS  coprime  moduli  [3  4  5]  (with  dynamic  range  Mms  =  60 )  is  chosen  to 
examine  the  idea  of  extending  the  unambiguous  range.  A  MATLAB  program  is  used  to 
plot  the  results  for  detecting  two  targets  that  lie  in  an  ambiguous  detection  range  bin  if  a 
single  modulus  is  used.  To  solve  the  unambiguous  range,  the  combined  residues  from  N 
sequences  are  used.  The  result  shows  that  the  paired  terms  from  both  targets  are  different, 
meaning  that  both  targets  have  different  corresponding  ranges.  Thus,  the  residue-Frank 
phase  code  sequences  can  be  used  to  extend  the  unambiguous  detection  range  from 
Ru  =  cT  /  2  =  ctbNc  /  2  to  Ru  RNS  =  ctbMms  /  2  where  c  is  the  speed  of  light  in  free  space. 

A  block  diagram  of  a  polyphase  CW  radar  system  that  uses  the  residue-Frank  phase  code 
is  shown  below  in  Figure  1. 
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Figure  1.  Block  diagram  of  a  CW  polyphase  radar  system  that  uses  a  residue -Frank 
phase  code. 


The  relationship  between  the  P4  phase  code  and  the  SNS  is  developed  based  on 
their  similar  parabolic  waveform  distribution.  A  new  equation  using  the  SNS  to  derive 
the  P4  phase  code  is  formulated.  This  new  phase  code  sequence  is  called  the  symmetrical 
residue-P4  phase  code.  The  PACF  and  the  PAF  show  that  the  symmetrical  residue-P4 
phase  code  formula  gives  the  same  results  as  using  the  formula  in  the  P4  phase  code. 
Then,  the  symmetrical  residue-P4  phase  code  sequences  using  a  set  of  N  =  3  SNS 
coprime  moduli  [7  8  9]  (with  dynamic  range  M SNS  =  37  )  are  examined  for  the  two  targets 

located  in  an  ambiguous  range  bin  by  using  MATLAB  programming.  Similar  to  the 
residue-Frank  phase  code  sequences,  the  result  shows  that  the  symmetrical  residue-P4 
phase  code  sequences  can  be  used  to  extend  the  unambiguous  detection  range  from 

Ru=cT 1 2  =  cthNc  /  2  to  RuSNS  =  ctbMSNS  /  2. 

Reconstruction  of  the  target’s  range  from  the  paired  RNS  and  SNS  phase  terms 

can  contain  inaccuracies  when  the  incoming  target  range  straddles  two  range  bins.  As  a 

result,  the  recovered  range  of  the  target  has  a  large  error.  To  prevent  the  target  range 

error,  the  RSNS  is  used.  The  paired  terms  from  all  N  sequences  in  the  RSNS  when 

considered  together  change  one  at  a  time  at  each  range  bin  subcode  transition  (integer 

Gray  code  property).  This  property  makes  the  RSNS  particularly  attractive  for  controlling 
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the  target  range  error.  In  this  thesis,  the  relationship  between  the  P4  phase  code  and  the 
RSNS  is  also  developed.  A  new  equation  using  the  RSNS  to  derive  the  P4  phase  code  is 
formulated.  This  new  phase  code  sequence  is  called  the  robust  symmetrical  residue-P4 
phase  code.  Within  each  P4  subcode,  the  phase  shifts  in  the  robust  symmetrical  residue- 
P4  phase  sequences  are  repeated  N  times  (equal  to  the  number  of  moduli),  resulting  in 
the  broader  mainlobe  for  the  PACF  and  the  reduction  of  signal  bandwidth  as  a  function  of 
N.  The  PACF  and  the  PAF  show  that  the  robust  symmetrical  residue-P4  phase  code 
formula  gives  similar  results  as  using  the  formula  for  the  P4  phase  code.  A  MATLAB 
program  is  used  to  examine  the  abilities  to  extend  the  unambiguous  range  and  to  control 
the  range  bin  errors  by  using  a  set  of  N  =  3  RSNS  coprime  moduli  [3  4  5]  (with  dynamic 
range  M RSNS  =  43 ).  The  results  verify  that  the  unambiguous  range  can  be  extended  from 

Ru=cT  12  =  ethNc  /  2  to  Ru  RSNS  =  ct'bM RSNS  /  2  where  t'b  is  the  subcode  period  in  the 
RSNS  phase  code,  and  the  range  bin  has  an  error  of  only  one  range  bin. 

In  summary,  all  three  NTT  are  investigated  to  show  that  they  can  be  used  to 
extend  the  unambiguous  range.  Also,  the  target  range  error  from  coding  can  be  prevented 
by  using  the  robust  symmetrical  residue -P4  phase  code  sequences.  The  use  of  these  new 
derived  phase  sequences  can  be  applied  to  CW  radar  systems  in  order  to  have  a  low 
probability  of  intercept  characteristic  given  that  the  set  of  moduli  and  the  radar  signal 
parameters  are  chosen  properly.  The  next  step  in  this  research  is  to  investigate  the 
performance  of  using  various  sets  of  coprime  moduli,  and  the  hardware  implementation 
and  testing  of  an  actual  system  and  is  being  proposed  for  future  work. 
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I.  INTRODUCTION 


A.  CONTINUOUS  WAVE  RADAR  SYSTEMS  USING  POLYPHASE 

MODULATION 

Continuous  wave  (CW)  radar  systems  have  a  superior  low  probability  of  intercept 
(LPI)  performance  over  pulse  train  radar  systems  due  to  their  low  average  power 
transmitted  (e.g.,  1  W)  and  their  use  of  pulse  compression  techniques.  Note  we  extend  the 
concept  of  pulse  compression  to  unpulsed  CW  waveforms  since  the  techniques  are 
similar  and  the  objectives  are  the  same.  Power  management  is  also  a  benefit  when  a  solid 
state  phased  array  antenna  is  used.  This  enables  the  target  return  signal-to-noise  ratio 
(SNR)  to  remain  at  a  constant  value  within  the  radar  receiver,  thus  lowering  the 
transmitted  power  further  as  the  range-to-target  decreases.  An  LPI  radar  system  is 
defined  as  a  radar  system  that  uses  a  special  emitted  waveform  intended  to  detect  targets 
as  well  as  to  prevent  a  non-cooperative  intercept  receiver  from  intercepting  and  detecting 
its  emission  [1].  On  the  other  hand,  pulsed  radar  systems  require  a  relatively  high  peak 
power  (e.g.,  10  kW)  to  be  transmitted  to  obtain  the  same  probability  of  target  detection. 
Consequently,  the  non-cooperative  intercept  range  is  significantly  longer  for  the  pulsed 
radars. 

Since  pure  CW  waveforms  cannot  resolve  the  target’s  range,  periodic  modulation 
techniques  are  used,  such  as  frequency  modulated  CW  (FMCW),  frequency-shift  keying 
(FSK),  noise  modulation,  phase-shift  keying  (PSK),  as  well  as  hybrids  of  these 
techniques.  The  first  step  to  designing  CW  radar  systems  that  use  periodic  modulations 
for  compression  is  to  decide  on  the  range  resolution  that  is  required.  This  in  turn  sets  the 
transmitted  bandwidth  of  the  waveform  for  the  above  techniques  (except  for  FSK  where 
the  range  resolution  is  dependent  on  the  duration  of  each  frequency).  Due  to  the  advances 
in  high-speed  processing  and  direct  digital  synthesis  modules  [2],  the  use  of  PSK 
techniques  in  CW  radar  is  highly  advantageous.  CW  radars  that  transmit  and  receive  PSK 
signals  can  result  in  LPI  radar  systems  with  small  range  resolution  cells  and  are  ideally 
suited  for  many  sensor  applications  for  situational  awareness,  including  multiple-input 
multiple-output  (MIMO)  configurations. 
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Binary  Barker  phase  coded  sequences  (1,  -1)  are  one  of  the  most  popular  PSK 
techniques  in  the  design  of  pulse  compression  radar  signals  due  to  the  aperiodic 
autocorrelation  coefficients  or  time  sidelobes  being  limited  to  1  /  Nc  (at  zero  Doppler) 

relative  to  a  mainlobe  level  of  one,  where  Nc  is  the  number  of  subcodes  used.  The 
longest  code  length  is  Nc=  13  (or  Nc=  169  for  a  compound  code)  [3],  These  codes  are  not 
used  in  CW  radar  systems  for  LPI  applications  since  they  are  very  sensitive  to  Doppler 
shifts  and  can  be  easily  detected  by  an  intercept  receiver  that  uses  frequency  doubling 
detection,  which  consists  of  multiplying  the  signal  by  itself  and  processing  the  result 
through  an  envelope  detector. 

Polyphase  sequences  are  also  of  finite  length  with  Nc  subcodes  and  consist  of 

discrete  time  complex  values  with  constant  amplitude  but  with  a  variable  phase. 
Polyphase  coding  refers  to  phase  modulation  of  the  CW  carrier  with  a  polyphase 
sequence  consisting  of  a  number  of  discrete  phase  states  with  each  phase  corresponding 
to  a  subcode.  The  number  of  subcodes  is  taken  from  an  alphabet  of  size  Nc  >  2,  which  is 
also  the  processing  gain  of  the  radar  excluding  any  post  detection  integration.  Increasing 
the  number  of  phase  values  in  the  sequence  allows  the  construction  of  longer  sequences, 
resulting  in  a  greater  processing  gain  or  compression  ratio  in  the  receiver  (or  equivalently 
a  larger  SNR).  The  code  period  is  T  =  Nctb  where  the  size  of  each  subcode  th  determines 
the  3  dB  bandwidth  of  the  wavefonn  B  =  \/tb  and  the  range  resolution  A R  =  cth  /  2.  The 
unambiguous  range  is  limited  by  the  number  of  subcodes  in  the  sequence  as  Ru  =  NCAR. 

With  a  CW  polyphase  waveform,  the  matched  filter  in  the  receiver  is  a  coherent, 
range-Doppler  matrix  correlation  processor  that  performs  a  cross  correlation  between  the 
received  signal  and  a  reference  signal  whose  envelope  is  the  complex  conjugate  of  Nr 
code  periods  of  the  transmitted  polyphase  signal.  The  cross  correlation  output  values  are 
added  together  to  reduce  the  ambiguity  sidelobe  levels.  When  the  polyphase  signals  are 
returned  and  have  an  impressed  Doppler  shift,  the  correlation  process  used  in  the  receiver 
compression  operation  is  not  perfect,  resulting  in  a  certain  amount  of  correlation  loss  due 
to  the  phase  shifts  across  the  code  period  T. 
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Two  of  the  most  important  PSK  codes  that  are  used  in  CW  radar  systems  include 
the  Frank  code  [4]  and  the  P4  code  [5],  The  Frank  code  is  a  polyphase  code  that  is 
derived  from  a  step  approximation  to  a  linear  FMCW  waveform  using  M  frequency 
steps  and  M  samples  per  frequency,  giving  the  number  of  subcodes  or  processing  gain  as 
Nc  =M2.  The  phase  steps  of  the  Frank  code  exhibit  a  saw-tooth  folding  waveform 
distribution.  The  peak  sidelobe  level  (PSL)  of  a  single  compressed  Frank  code  period  is 
PSL  =  201og10(l/(M;r)).  The  P4  code  consists  of  discrete  phases  of  a  linear  chirp 

waveform  taken  at  specific  time  intervals.  It  is  derived  by  converting  a  linear  FMCW  to 
baseband  using  a  local  oscillator  on  one  end  of  the  frequency  sweep  and  sampling  the 
inphase  and  quadrature  video  at  the  Nyquist  rate.  The  P4  phase  steps  exhibit  a 
symmetrical  parabolic  distribution.  The  PSL  of  a  single  P4  compressed  code  period  is 

PSL  =  201og10  yj2/  (Nc7r2). 

Although  polyphase  CW  modulations  can  serve  as  a  transmission  waveform  for  a 
LPI  radar  system,  there  are  several  limiting  factors.  Most  significantly,  the  unambiguous 
detection  range  of  the  waveform  is  limited  by  the  number  of  subcodes  (code  length)  as 
Ru  =  NCAR.  By  increasing  the  number  of  subcodes  to  extend  the  unambiguous  range,  a 
larger  range-Doppler  correlation  matrix  processor  in  the  radar  receiver  is  needed  to 
provide  the  code  compression.  To  lower  the  peak  side  lobes,  Nr  copies  of  the  phase  code 
are  used  to  compress  the  N  returned  code  periods  from  the  target.  Consequently,  the 

increased  code  length  or  processing  gain  requires  a  larger  code  compression  time  as  well 
as  an  increase  in  the  bulk  memory  requirements  in  the  receiver.  Further,  a  significant 
amount  of  correlation  loss  is  incurred  if  the  total  number  of  code  periods  returned  from 
the  target  N  is  less  than  the  number  of  code  periods  Nr  used  in  the  correlation 

processor,  as  would  be  the  case  due  to  a  limited  time-on-target. 

In  this  thesis,  a  novel  relationship  between  the  Frank  code  and  the  residue  number 
system  (RNS)  is  developed  in  order  to  significantly  extend  the  unambiguous  range 
beyond  a  single  code  period  T  =  Nctb  while  not  having  to  increase  the  number  of 

subcodes.  In  addition,  due  to  the  symmetrical  distribution  of  the  P4  polyphase  code,  the 
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unambiguous  range  is  extended  by  developing  a  phase  relationship  based  on  symmetrical 
number  systems;  both  the  symmetrical  number  system  (SNS)  and  the  robust  symmetrical 
number  system  (RSNS)  are  investigated.  The  phase  relationships  are  constructed  using 
these  number  theoretic  transforms  to  extend  the  unambiguous  range  of  the  radar  system 
by  associating  each  phase  subcode  within  the  code  period  T  with  an  integer  sequence 
within  a  coprime  modulus  from  the  set  .  The  paired  integer  terms  from 

each  polyphase  sequence  are  unambiguous  within  the  transform’s  dynamic  range.  The 
dynamic  range  is  defined  as  the  greatest  length  of  combined  phase  sequences  that  contain 
no  ambiguities  or  repeated  paired  terms.  The  dynamic  range  for  the  RNS  is  denoted 

Mms  (Frank  code).  The  dynamic  range  for  the  SNS  is  MSNS  (P4  code).  The  dynamic 
range  for  the  RSNS  is  MRSNS  (P4  code). 

For  each  case,  the  length  of  the  dynamic  range  contains  multiple  code  periods.  In 
the  radar  receiver,  the  integer  values  within  each  modulus  sequence  correspond  to  the 
phase  values  from  a  phase  detector  [6],  The  RNS  phase  sequence  contains  multiple  Frank 
codes  and  extends  the  unambiguous  target  detection  range  from  Ru=cT  /  2  =  cthNc  /  2  to 
R„  rns  =  ctb( n /«, )  /  2  =  ctbMms  /  2.  The  SNS  phase  sequence  contains  multiple  P4  codes 
and  extends  the  unambiguous  target  detection  range  beyond  a  single  code  period  to 
Ru  SNS  =  cthMSNS  /  2.  It  is  noted  that  the  SNS  also  has  a  direct  relationship  with  the 

discrete  Fourier  transform  (DFT)  and  has  been  used  to  resolve  frequency  ambiguities  in 
undersampling  digital  receivers  [7].  Similarly,  the  unambiguous  range  using  the  RSNS  is 

extended  to  Ru  RSNS  =  ct'hMRSNS  /  2  where  t'b  is  the  subcode  period  in  the  RSNS  phase 

code. 

Reconstruction  of  the  target’s  range  from  the  paired  RNS  and  SNS  phase  terms 
can  contain  inaccuracies  when  the  incoming  target  range  straddles  two  range  bins  or 
subcodes.  For  example,  if  the  target  is  straddling  a  range  bin  and  a  constant  false  alarm 
rate  (CFAR)  processor  detects  the  target  in  the  wrong  subcode  within  a  modulus 
sequence,  the  recovered  range  of  the  target  will  have  a  large  error.  This  is  shown  to  be  the 
case  for  both  the  RNS  and  the  SNS.  Additional  signal  processing  can  be  used  to  eliminate 
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these  errors  but  adds  additional  complexity  to  the  detection  process.  A  more  efficient 
method  to  eliminate  the  target  range  bin  errors  uses  the  RSNS.  The  RSNS  is  also  a 
modular  scheme;  however,  the  paired  terms  from  all  N  sequences,  when  considered 
together,  change  one  at  a  time  at  each  range  bin  subcode  transition  (integer  Gray  code 
property).  Although  the  dynamic  range  of  the  RSNS  is  less  than  the  SNS 
(Mrsns  <  Msns  ),  the  RSNS  integer  Gray  code  property  makes  it  particularly  attractive  for 

controlling  the  target  range  error.  In  summary,  all  three  number  theoretic  transforms  are 
investigated  to  determine  the  feasibility  of  extending  the  unambiguous  range  of  the  CW 
polyphase  emitter. 

B.  PRINCIPAL  CONTRIBUTIONS 

The  first  step  in  this  thesis  was  to  understand  the  principles  of  CW  phase 
modulation  radar  systems.  This  included  the  complex  envelope  representation  of  the 
transmitted  signal  and  the  related  modulation  parameters  such  as  the  code  period,  the 
subcode  period,  the  number  of  subcodes  within  the  code  period,  and  the  number  of  cycles 
of  the  carrier  frequency  per  subcode.  Their  relationship  to  the  signal  bandwidth,  the  range 
bin  size,  the  unambiguous  detection  range  and  the  processing  gain  were  also  studied. 
Upon  reception,  the  CW  signal  is  compressed,  and  several  code  compression  techniques 
were  studied  and  modeled.  To  understand  the  receiver  performance  when  a  particular 
CW  phase  code  is  used,  the  periodic  autocorrelation  function  and  the  periodic  ambiguity 
function  were  examined.  A  MATLAB  code  was  written  to  generate  the  CW  phase  coded 
signals  and  to  call  the  periodic  ambiguity  function  subroutines. 

The  next  step  in  the  thesis  was  to  study  the  Frank  polyphase  code  formulation  and 
the  phase  distribution  within  a  code  period.  Periodic  ambiguity  analysis  was  also 
performed  to  quantify  the  sidelobe  levels  within  the  receiver’s  ambiguity  space  (range 
offset  -  Doppler  offset).  Next  the  RNS  was  studied  as  a  means  to  generate  the  Frank  code 
since  they  both  have  sawtooth  folding  waveform  representations.  By  using  the  RNS  to 
create  the  Frank  polyphase  modulation,  the  unambiguous  range  of  the  code  can  be 
extended  significantly  beyond  a  single  code  period.  A  new  expression  to  generate  the 
Frank  polyphase  code  was  developed  as  a  function  of  the  coprime  RNS  moduli  chosen. 
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By  using  two  or  more  Frank  code  sequences  that  are  relatively  prime,  the  detected  target 
residue  from  each  sequence  is  combined  within  a  system  of  congruences  which  can  be 
solved  in  a  straightforward  manner  using  the  Chinese  Remainder  Theorem  (CRT).  The 
solution  is  the  target’s  unambiguous  range.  A  MATLAB  model  was  developed  to 
simulate  the  new  Frank  polyphase  modulation  in  order  to  verify  the  expression  and 
determine  the  differences  in  phase  values  (if  any)  that  might  occur.  Next  a  new  CW  radar 
signal  processing  algorithm  was  developed  to  transmit  and  receive  the  coprime  Frank 
sequences  to  demonstrate  the  extension  of  the  unambiguous  range.  The  algorithm  also 
addressed  multiple  targets  within  the  field  of  view.  Several  examples  were  given  to 
demonstrate  the  concept.  To  examine  the  practical  considerations,  a  simulation  of  the 
transmit  power  versus  required  compression  output  SNR  was  performed  in  MATLAB  for 
several  values  of  coprime  moduli. 

Next,  the  P4  polyphase  code  was  examined  the  relationships  to  the  SNS  and  the 
RSNS.  A  new  formulation  of  the  P4  polyphase  code  as  a  function  of  the  SNS  and  RSNS 
coprime  moduli  was  developed.  Analysis  similar  to  that  for  the  Frank  code  described 
above  was  performed  to  verify  the  new  P4  expression.  Finally,  a  summary  of  the  different 
relationships  that  were  developed  is  given  and  a  comparison  was  perfonned. 

C.  THESIS  OUTLINE 

In  Chapter  II,  the  fundamentals  of  radar  systems  that  use  polyphase  modulation 
are  reviewed.  In  Chapter  III,  the  Frank  code  and  RNS  are  examined  to  extend  the 
unambiguous  range.  In  Chapter  IV,  the  relationships  between  the  P4  and  the  SNS  are 
developed.  In  Chapter  V,  the  relationships  between  the  P4  and  the  RSNS  are  investigated 
to  reduce  the  possibility  of  a  range  bin  encoding  error.  Concluding  remarks  are  offered  in 
Chapter  VI. 
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II.  PRINCIPLES  OF  PHASE  MODULATION  RADAR 


The  idea  of  CW  radar  using  phase  modulation  is  reviewed  in  this  chapter.  First, 
the  fundamentals  and  use  of  phase  modulation  in  CW  radar  systems  are  explained.  Next, 
the  concepts  of  signal  analysis  and  the  technique  used  in  the  receiver’s  signal  processing 
are  presented.  Also,  the  system  equations  for  a  phase  modulation  in  CW  radar  are  given 
in  order  to  examine  the  signal’s  properties.  These  equations  will  be  used  to  determine  the 
new  phase  modulation  code’s  performance  in  the  following  chapters. 

A.  CW  PHASE  MODULATION  CONCEPTS 

In  Figure  1,  the  basic  CW  radar  geometry  is  illustrated.  First,  the  CW  radar 
transmits  the  phase  coded  signal  s(t )  to  the  target  through  the  transmitting  antenna. 
When  the  signal  propagates  to  the  target  at  a  range  R,  it  reflects  back  to  the  CW  radar 
with  different  characteristics.  For  example,  the  signal’s  amplitude  is  less  and  the 
frequency  is  shifted  in  time  in  accordance  with  the  relative  movement  between  platforms 
along  the  line-of-sight  (V cos#  in  Figure  1)  or  Doppler  shift  v.  The  return  signal  x(t)  is 
then  captured  by  the  receiving  antenna.  In  the  radar  system’s  receiver  and  digital  signal 
processor,  the  compression  processing  output  is  used  to  detect  and  extract  the  target’s 
characteristics. 

In  phase  coded  CW  radar  systems,  the  phase  shifting  operation  is  perfonned  in  the 
radar’s  transmitter,  with  the  timing  information  generated  from  the  receiver-exciter.  The 
transmitted  complex  signal  can  be  written  as  [1] 

s(t)  =  Aej(2xf't+A(t))  (1) 

where  A  is  the  signal  amplitude,  fc  is  the  carrier  frequency,  and  </>k(t)  is  the  time 
dependent  phase  modulation  code.  The  inphase  I  and  quadrature  Q  representations  of 
the  complex  signal  from  the  transmitter  can  be  represented  as 

I  (0  =  A  cos(2;r  fct  +  <f>k  (t))  (2) 


and 


Q(t)  =  Asin(27rfct  +  </>k(t)). 
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Moving 


Target 


CW  phase 
modulation 
radar 


Phase  code 


Transmitter 


Code 

compression 


Figure  1 .  Geometry  of  CW  radar  using  phase  modulation  codes  to  sample  a  moving 

target. 


For  one  code  period  T,  the  CW  signal  is  shifted  in  phase  every  subcode  with  each 
code  period  consisting  of  Nc  subcodes.  Each  subcode  with  phase  </>k  has  a  duration  of  tb 
which  is  the  subcode  period.  For  a  specific  code  sequence,  the  code  period  is 


T  =  Nctb. 


The  transmitted  signal  can  be  expressed  as 


=  YJuk[t-{k-\)tb] 


k= 1 


for  0  <  t  <  T  and  zero  elsewhere.  The  complex  envelope  uk  is 


(4) 

(5) 

(6) 


If  cpp  is  the  number  of  cycles  of  the  carrier  frequency  per  subcode,  the  bandwidth  B  of 
the  transmitted  signal  is 

fa  1 


B  = 


cpp  tb 


(V) 
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To  understand  the  concepts  of  phase  modulation,  a  binary  Barker  phase  code  is 
examined  [3].  In  Figure  2,  an  example  of  one  code  period  with  Nc=l  is  shown.  Here, 

the  carrier  frequency  is  fc  =  1  kHz  ,  the  subcode  period  th  =  1  ms ,  and  from  (7)  cpp  =  1 . 
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Figure  2.  Phase  coded  signal  with  the  Barker  code  Nc=l ,  A  =  1 ,  fc  =  1  kHz  , 

tb  =  1  ms. 


The  received  return  signal  x(t)  from  the  target  can  be  written  as 

x(t )  =  7^{2jr(/'+vX'-r)+^(,-r)}  (8) 

where  //  is  the  amplitude  attenuation  coefficient  of  the  received  signal  due  to  target 
scattering  and  propagation  attenuation.  The  two-way  roundtrip  delay  r  for  a  monostatic 
radar  system  is  t  —  2R/c,  where  c  is  the  speed  of  light  in  free  space.  The  Doppler 
frequency  v  due  to  the  target  motion  is  v  =  2V  cos(0)  /  A  where  A  is  the  signal 
wavelength  and  V  is  the  target  velocity. 

B.  CODE  COMPRESSION  OF  RECEIVED  SIGNAL 

Code  compression  (CC)  is  a  signal  processing  technique  designed  mainly  to 
minimize  the  range  resolution  cell  and  increase  the  SNR  of  the  returned  signal.  The  idea 
of  CC  is  similar  to  that  in  signal  correlation  in  pulse  train  radar.  As  shown  in  Figure  3  a, 
the  CW  radar  repeats  the  phase  sequence  every  code  period  T.  A  similar  pulse  train 
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signal  would  also  have  a  pulse  repetition  interval  PRI  =  T.  Referring  to  Figure  3b,  if  the 
pulse  width  rpw  of  the  pulse  train  signal  is  equal  to  the  width  of  the  subcode  period  tb  in 

the  CW  phase  coded  signal,  the  output’s  mainlobe  from  both  correlator  and  compressor 
would  have  the  same  duration  ( 2 r  =  2 th  ).  These  CC  peaks  represent  the  targets  if  the 

peaks  are  higher  than  the  threshold  desired  in  a  constant  false  alarm  rate  (CFAR)  process. 
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Figure  3.  The  relationship  between  uncompressed  pulse  train  and  phase-coded  CW 
signal  for  (a)  Transmitted  signal,  and  (b)  The  output  of  the  signal  processed  in  the 

receiver. 


The  phase-coded  signal  is  transmitted  to  detect  the  target  at  a  range  R.  When 
reaching  the  target,  the  reflected  signal  returns  to  the  radar.  At  the  radar  receiver,  the 
phase  coded  return  signal  is  correlated  with  the  reference  phase  code.  This  is  achieved  by 
convolving  the  incoming  signal’s  phase  code  with  a  conjugated  and  time-reversed 
version  of  the  transmitted  signaFs  phase  code.  This  operation  can  be  done  either  in  the 
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digital  or  analog  domain.  Many  digital  signal  processing  architectures  have  been 
proposed  for  fast  correlation  using  fast  Fourier  transforms.  If  no  compression  loss  is 
present,  the  correlated  output  will  give  a  very  narrow  pulse  width  at  the  compression 
output  for  the  detected  target.  The  correlation  output  is  then  sent  to  a  target  detection 
process  such  as  a  constant  false  alarm  rate  processor.  If  there  is  any  alteration  of  the 
return  signal’s  phase  (for  example  from  a  Doppler  shift  or  a  less  than  adequate  bandwidth 
in  the  receiver  B  <1  !tb  ),  a  certain  degree  of  correlation  loss  occurs. 

The  received  waveform  from  the  target  is  digitized  and  correlated  using  a  phase 
code  compressor.  The  range-Doppler  correlation  matrix  can  contain  a  cascade  of  Nr  sets 

of  Nc  reference  coefficients  to  improve  the  sidelobe  structure.  In  Figure  4,  the  range- 
Doppler  correlation  matrix  receiver  for  zero  Doppler  offset  is  illustrated.  Here,  the 
received  signal  consists  of  N  code  periods,  each  of  which  has  Nc  subcodes.  The  return 

signal  is  first  processed  by  a  filter  matched  to  a  rectangular  subcode  of  length  tb , 
followed  by  a  phase  detector  that  sends  forward  the  phase  values.  The  detected  output 
signal  is  then  sent  through  a  tapped  delay  line  where  each  delay  D  is  tb.  As  the  signal 
progresses  through  the  tapped  delay  line,  it  is  multiplied  by  the  reference  signal.  At  each 
step,  the  multiplication  for  each  delay  is  summed  separately  for  each  of  the  Nr  reference 

code  periods.  To  reduce  the  sidelobe  structure,  a  weighting  function  C(.  can  also  be  added 
[8].  To  compute  the  entire  range-Doppler  ambiguity  function,  the  correlation  outputs 
from  each  subcode  are  multiplied  by  qk  =  e’2akAvtb  for  v  =  A  v  where  k  ranges  from  0  to 
NrNc  - 1 .  For  v  =  2Av,  the  output  of  the  previous  v  =  A  v  operation  are  again  multiplied 
by  qk  =  ej2rrkAv‘h  and  so  on  for  the  entire  ambiguity  space  [9]. 


11 


Correlation  Output 


Figure  4.  Correlation  receiver  matched  to  Nr  periods  of  a  transmitted  polyphase  code 

for  v  =  0  (After  [9]). 

In  Figure  5,  an  illustration  of  the  code  compression  for  CW  binary  Barker  code 
with  Nc  =  7  is  shown.  The  radar  receives  the  returned  Barker  sequence  with  the  number 

of  code  periods  N  =3.  The  number  of  reference  code  periods  used  in  the  receiver  is 
Nr  =  1 .  The  receiver  starts  compressing  the  first  returned  code  Tx  at  r  =  0,  giving  the 
correlation  output  peak  at  a  normalized  delay  of  r  /  tb  =  0  representing  the  detected 
target.  At  each  step  in  the  compression,  the  returned  signal  is  shifted  one  subcode  period 
r  =  tb.  When  the  returned  code  is  shifted  by  one  code  period  T,  the  correlation  output 

gives  another  peak  value  due  to  the  second  code  T2  being  compressed  at  a  normalized 
delay  of  r  I  tb  =  7.  The  last  peak  from  the  third  code  T3  compressed  is  also  shown  at  a 
normalized  delay  offset  of  r  /  tb  =  14.  That  is,  the  peak  value  for  the  detected  target 
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repeats  every  code  period  T.  Generally,  the  width  of  the  PACF  mainlobe  is  2 tb.  To  be 
able  to  distinguish  between  two  targets  close  to  one  another,  the  pulses  have  to  be 
separated  by  at  least  r  =  tb  period  (which  is  defined  as  the  range  resolution  cell). 


Reference  code 


T, 

JL 


Returned 

code 


T, 

JL 


r  -0;  PACF-0  xl Wlx  lHlx  l)+<  -lx  -l)+(  -lx-l)-t(lx  l)+(  -lx  -1).  7 


r=l,  PACF=(lxl)4<lxl)-Klx-l)+(-lx-l)4(-lxl)-Klx-l)4(-lxl)=-l 


*  r  =  X  PACF=(lxl)-l(lx-l)-Klx-l)4<-lxl)4(-lx-l)4(lxl)-K-lxl)=-l 
r  =  3,  PACF-(lx-l)-t<lx-l)-t<lxl)+(:-lx-l)'K-lxO'Klxl)'K-lxl)=-l 
^  ^  r  =  4,  PACF=(lx-l)4(lxl)+(lx-l)+(-lxl)+(-lxl)4<lxl)-K-lx-l)--l 


* 

IF 

* 

* 

* 

IF 

* 

+ 

+  '  |  r  =  5,  PACF=(lxl)-Klx-l)4<lxl)+(-lxl)-K-lxl)4<lx-l)4<-lx-l)=- 

+ 

F 

♦ 

+ 

* 

T7 

♦ 

+ 

+ 

|  ^r-3  r=  6,  PA2F=(lx-l)+<lxl)+(lxl)+(-lxl)+(-lx-l)+(lx-l)+(-lxl)=-l 

^  1  =  7,  PACT=<lxl)*<lxl)t<lxl)4<-lx-l)*<-lx-lXlxlK-lx-l)=7 

* 

* 

♦ 

* 

♦ 

* 

* 

* 

* 

♦ 

♦ 

♦ 

IE 

♦ 

♦ 

♦ 

V-1  r  =  8,  PACF=(lxl)-t<lxl)-Klx-l)-K-lx-l)-K-lxl)-t<lx-l)-K-lxl)=-l 

Figure  5.  Example  of  code  compression  in  CW  Barker  code  Nc=l  and  N  =3. 


In  summary,  the  unambiguous  range  of  a  specific  CW  phase-coded  sequence  is 

*-=?=£ N~r  (9> 

and  the  range  resolution  for  one  transmitted  polyphase  signal  is 


A R  =  «L. 
2 


(10) 


Another  important  parameter  in  pulse  compression  is  the  pulse  compression  ratio 
(PCR)  or  time -bandwidth  product.  This  is  also  referred  to  as  the  radar’s  processing  gain 
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(PG)  and  is  the  preferred  term  in  this  thesis.  In  general,  the  higher  the  PG,  the  higher  the 
SNR  improvement  in  the  signal  processing  step.  The  radar  signal  processing  input  SNR  is 
SNRm  and  the  radar  signal  processing  output  SNR  is  SNRRo( see  Figure  1).  They  are 
related  to  the  processing  gain  of  the  signal  processor  by 

pg  =  ^r!L  (H) 

SNRRi 

For  a  polyphase  code,  the  PG  is  equal  to  the  time-bandwidth  product  or  the  number  of 
subcodes  as 

PG  =  T(l/tb)  =  (Nctb)/tb=Ne.  (12) 

In  the  Barker  code  Nc  =  7  example,  the  processing  gain  is  7.  For  a  carrier 
frequency  fc,  narrowing  the  subcode  period  th  enables  the  transmitted  signal  to  be 
spread  over  a  larger  bandwidth.  Also,  increasing  the  number  of  subcodes  Nc,  the 

correlation  output  sidelobes  can  be  decreased  significantly  while  providing  a  larger  PG. 
This  comes  at  the  expense  of  a  larger  range-Doppler  correlation  matrix  processor  in  the 
receiver.  In  addition,  a  longer  compression  time  and  an  increase  in  the  size  of  the  bulk 
memory  is  required. 

The  signal  processing  function  used  for  the  compression  is  the  autocorrelation 
function  (ACF).  For  CW  radars,  the  compression  is  accomplished  with  a  periodic 
autocorrelation  function  (PACF).  Also,  the  periodic  ambiguity  function  (PAF)  is  used  to 
show  the  delay-offset  versus  frequency-offset  in  the  receiver  and  shows  the  resulting 
sidelobe  structure.  This  is  especially  useful  for  monitoring  the  ambiguity  space  for 
secondary  targets  that  might  hide  in  the  sidelobe  structure.  Note  that  the  zero  Doppler  cut 
of  the  PAF  is  the  PACF  [1], 

C.  PERIODIC  AUTOCORRELATION  FUNCTION 

It  is  important  that  in  order  to  achieve  high  range  resolution  for  a  polyphase  code 
sequence,  the  coded  signal  has  to  have  a  noise-like  PACF  property,  close  to  an  impulse 
function.  The  ACF  is  used  for  signals  of  finite  energy.  For  a  CW  polyphase  code 
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sequence  with  the  number  of  phase  codes  Nc ,  subcode  duration  tb ,  and  a  code  period 
T  =  Nctb ,  a  periodic  complex  envelope  u(t )  is  given  by  [1] 

u(t)  -  u{l  +  nT )  (13) 

for  n  -  0,±l,±2,±3,....The  values  of  the  PACF  as  a  function  of  the  delay  r  (which  are 
the  multiples  of  tb )  are  given  by 


R(rth)  =  ^-YJu{n)u(n  +  r).  (14) 


Ideally,  a  perfect  periodic  autocorrelation  function  (zero-sidelobe)  is  desired  where 
R(rth)  =  1  when  r  =  0  (mod  A  ),  and  R{rtb)  =  0  when  r  ^  0  (mod  Nc).  Since  the  CW 
signal  is  continuous,  a  perfect  periodic  autocorrelation  function  can  be  achieved. 


In  Figure  6a,  the  ACF  output  for  a  Nc  =  7  Barker  code  is  plotted  (in  dB)  with  no 
noise  using  fc  =  1  kHz ,  fs  =  7  kHz  ,  cpp  =  1  and  Nr  =  1  using  the  MATLAB  m-fde 
(amfbn7.m)  [11].  The  highest  peak  mainlobe  is  at  the  zero  time  delay.  The  width  of  the 
mainlobe  is  equal  to  the  length  of  the  subcode  tb.  For  the  Barker  code,  it  is  found  that  the 

PSL  is  equal  to  -201ogl0(Ac)  =  -17  dB.  The  PACF  in  Figure  6b  shows  that  the  peak  is 
repeated  every  code  period  T  =  Nctb  at  the  value  of  7.  Also  note  from  (14)  that  the  CW 
binary  Barker  code  is  not  a  perfect  code. 


(a) 


(b) 


Figure  6. 


Normalized  correlation  output  of  Barker  code  sequence  Nc  =  7  and  N r  =  1 
for  (a)  ACF  and  (b)  PACF. 
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Today,  many  studies  on  PSK  show  that  the  PSL  of  the  phase-coded  sequences  can 
be  reduced  by  choosing  a  proper  phase  coding.  Moreover,  there  are  CW  polyphase  codes 
able  to  achieve  a  zero  sidelobe  (perfect  codes). 


D.  PERIODIC  AMBIGUITY  FUNCTION 

In  a  realistic  environment,  the  returned  signal  will  be  shifted  in  frequency  if  the 
target  is  moving  relative  to  the  radar’s  antenna  beam  axis  (Doppler  effect).  The  PAF  is  an 
important  tool  to  assess  the  response  of  a  correlation  receiver.  In  pulsed  radar  analysis, 
the  ambiguity  function  is  a  two-dimensional  function  of  time  delay  offset  z  and 

Doppler  frequency  offset  v  showing  the  correlation  between  the  returned  signal  and  the 
receiver  matched  filter.  For  the  finite  duration  signal,  the  ambiguity  function  (AF)  is 
determined  by  the  properties  of  the  pulse  and  the  matched  filter.  If  u{t)  is  the  complex 
envelope  of  both  the  transmitted  signal  and  received  signal,  the  AF  is  given  by  [10] 


|j(r,v)|=  J  u(t)u* (t -r)eJ2ml dt 


(15) 


where  *  denotes  the  complex  conjugate.  A  target  further  from  the  radar  than  the  reference 
(r  =  0)  position  will  correspond  to  positive  r.  A  positive  v  implies  a  target  moving 
toward  the  radar. 


A  more  concise  way  of  representing  the  ambiguity  function  consists  of  examining 
the  one-dimensional  zero-delay  and  zero-Doppler  cuts;  that  is,  j(0,u)  and  /(r.0) , 
respectively.  The  matched  filter  output  as  a  function  of  a  time  (the  signal  one  would 
observe  in  a  radar  system)  is  a  delay  cut,  with  constant  frequency  given  by  the  target’s 
Doppler  shift  x(r,v).  Ideally,  the  ambiguity  diagram  should  show  the  diagonal  ridge 
centered  at  the  origin  and  zero  elsewhere,  no  ambiguity. 

The  PAF,  introduced  by  Levanon  and  Freedman  [11],  describes  the  response  of  a 
correlation  receiver  to  a  CW  signal  modulated  by  a  periodic  waveform  with  period 
T,  when  the  reference  signal  is  constructed  from  an  integral  number  Nr  of  periods  of  the 

transmitted  signal  (coherent  processor  length  NrT ).  The  target  illumination  time  (dwell 
time)  NT  must  be  longer  than  NT.  As  long  as  the  delay  z  is  shorter  than  the 
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difference  between  the  dwell  time  and  the  length  of  the  reference  signal 
0  <  t  <  (N  - Nr)T,  the  illumination  time  can  be  considered  infinitely  long  and  the 

receiver  response  can  be  described  by  the  PAF  given  as  [8] 


NrTr 

J  u(t  -  t)u* (t)ei2nvt dt 
o 


(16) 


Figure  7.  PAF  of  Barker  code  Nc  =  7  with  (a)  Nr  =  1  and  (b)  Nr  =  4. 
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In  Figure  7a,  the  PAF  of  the  Barker  code  Nc  =  7  with  Nr  =  1  is  plotted  using 
MATLAB  m-file  (amfbn7.m)  from  [10].  Also,  the  PAF  of  the  Barker  code  Nc=l  with 
Nr=  4  is  plotted  in  Figure  7b.  By  increasing  Nr,  the  Doppler  sidelobes  have  been 
significantly  lowered. 

In  the  next  chapter,  the  PACF  and  PAF  are  used  to  characterize  the  CW  Frank 
code.  The  residue  number  system  is  introduced  and  a  relationship  to  the  Frank  code  is 
developed. 

Since  there  are  many  variables  used  throughout  this  thesis,  all  of  important 
variables  are  listed  in  the  Appendix  to  provide  unambiguous  meanings. 
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III.  FRANK  POLYPHASE  MODULATION  AND  THE  RESIDUE 

NUMBER  SYSTEM 

In  this  chapter,  the  Frank  phase  code  and  the  residue  number  system  (RNS)  are 
introduced.  Their  relationship  is  presented. 

A.  FRANK  PHASE  CODE 

The  Frank  code  is  derived  from  a  step  approximation  to  a  linear  frequency 
modulation  waveform  using  M  frequency  steps  and  M  samples  per  frequency  [4],  If  i 
is  the  sample  number  in  a  given  frequency  and  j  is  the  frequency  number,  the  phase  of 
the  z'th  sample  of  the  /th  frequency  is  given  by 

A  =<!>,, i  “^jO-OO'-O  (17) 

M 

where  k  is  the  time  index  for  a  single  phase  change  per  subcode  and  it  ranges  from  1  to 
M2,  and  = 

The  phase  values  and  the  phase  values  modulo  2 n  of  the  Frank  code  for 
M  =  8,  Nc  =  64  are  shown  in  Figures  8  and  9,  respectively.  The  carrier  frequency  is 

fc  =  1  kHz,  the  sampling  frequency  is  fs  =  7  kHz,  and  the  number  of  carrier  cycles  per 
subcode  is  cpp  =  1.  The  ACF  and  the  PACF  are  shown  in  Figure  10  where  Nr  is  the 
number  of  reference  code  periods  in  the  receiver.  The  delay  axis  is  normalized  by  the 
subcode  period  tb.  Note  from  Equation  (14)  and  Figure  10b  that  the  Frank  code  is  a 

perfect  code.  The  PSL  is  PSL  =  201og(l/(M;r))  or  -28  dB.  The  PAF  in  Figure  11 

shows  the  delay-offset  and  Doppler-offset  sidelobes.  Note  that  the  Doppler  tolerance  is 
reflected  in  the  way  the  sidelobes  are  arranged  about  the  mainlobe.  The  study  on  Doppler 
tolerance  is  presented  in  [5]  and  [12]. 
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40 


Figure  8.  Frank  code  modulation  for  M  =  8 ,NC  =  64  (From  [1]). 


Figure  9.  Signal  phase  (radians)  modulo  2 n  versus  k  -index  for  phase  change  of  the 

Frank  code  with  M  =  8,  Nc  =64  (From  [1]). 
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Figure  10.  ACF  (a)  and  PACF  (b)  for  the  Frank  code  with  M  =  8,  cpp  =  1,  and  N r  =  1 

(From  [1]). 


Figure  1 1 .  PAF  for  the  Frank  code  modulation  with  M  =  8,  cpp  =  1,  and  Nr  =  1 

(From  [1]). 
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B. 


THE  RESIDUE  NUMBER  SYSTEM 


The  RNS  was  introduced  by  Szabo  [13].  It  has  been  established  as  an  important 
tool  in  parallel  processing  applications  and  high  speed  computations  [  14]— [  1 6] .  The  RNS 
can  serve  as  a  source  for  extending  the  unambiguous  range  by  decomposing  the 
unambiguous  range  into  a  number  of  parallel  sub-ranges  (moduli)  that  have  a  smaller 
number  of  subcodes.  Each  sub-range  for  a  different  coprime  modulus  requires  only  a 
code  period  in  accordance  with  that  modulus.  A  much  larger  unambiguous  range  is 
achieved  after  the  results  of  these  smaller  sub-range  operations  are  recombined.  Thus,  by 
developing  a  relationship  between  the  RNS,  the  polyphase  distributions  within  the  Frank 
code  period,  and  the  target  detection  processing,  the  extension  of  the  unambiguous  range 
beyond  a  single  code  period  is  feasible. 

The  RNS  residue  distribution  is  based  on  a  saw-tooth  folding  waveform,  where 
the  discrete  residue  values  rise  gradually  and  fall  quickly  as  the  input  value  rises 
gradually  (i.e.,  in  mod  5  the  succession  of  discrete  values  for  an  increasing  input  is  0,  1, 
2,  3,  4,  0,  1,  2,  3,  ...).  A  single  RNS  sequence  can  be  generated  as 

a  =  Rm  (mod  m)  ( 1 8) 

where  0  <  Rm  <  m,  m  is  the  modulus,  and  Rm  is  the  residue  of  a  modulo  m.  The  integer 
values  Rm  within  each  modulus  sequence  are 

Rm  =  [0,l,...,m  — 1,]  (19) 

and  are  representative  of  the  saw-tooth  folding  waveform.  The  sequence  is  repeated  in 
both  directions  forming  a  periodic  sequence  with  the  period 

Prns  =  m-  (2°) 

Given  a  set  of  RNS  coprime  moduli,  a  complete  residue  number  system  is 
formed.  A  set  of  integers  y0,yl,...,ym_l (mod m)  form  a  complete  residue  system  if 

y  =  g(mod  m)  for  g  =  0,1,...,  m  -1 .  The  dynamic  range  Mms  of  the  RNS  is  the  number 
of  paired  terms  from  the  sequence  of  values  that  contain  no  ambiguities  and  is 
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(21) 


N 

^  RNS  —  I~ [  mi 
1=1 

where  mi  is  the  zth  coprime  modulus  and  N  >  2  is  the  number  of  moduli,  1  <  i  <  N.  In 
Figure  12,  an  RNS  example  for  N  =  2  moduli,  mx  =4  and  m2  =5,  is  shown.  Note  the 
stair  step  relationship  of  the  residues. 


inputs 


Figure  12.  The  RNS  residue  folding  waveforms  for  ml  =  4  and  m2  =  5 . 


From  Equation  (21),  for  ml  =  4  and  m2  =  5,  the  dynamic  range  is  M RNS  =  20.  In 
Table  1,  the  two  sequences  are  shown  along  with  the  input  a.  Starting  at  a  =  0,  the  total 
number  of  distinct  paired  terms  is  20,  0  <  a  <  19.  The  index  value  /  =  a  +  1  in  Table  1  is 
shown  to  count  the  integer  number  of  distinct  paired  terms  within  the  dynamic  range. 


Table  1 .  Finding  the  RNS  dynamic  range  for  ml  =  4  and  m2=  5 . 
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c. 


RESIDUE-FRANK  PHASE  CODE 


The  RNS  residues  Rm  k  can  be  folded  into  the  Frank  phase  code  sequence  by 
associating  Rm  k  with  the  Frank  subcode  phases.  The  phase  index  k  starts  from  1  to  Nc, 
where  Nc  is  the  number  of  subcodes  within  one  code  period.  The  residue -Frank  phase 
code  sequence  from  modulo  m  is  given  by 


,  2  n 

Ym,k  ~ 

m 


k- 1 


m 


R 


(22) 


where  k  =  1,2,3, ...,/m2  is  the  phase  index  for  each  phase  change,  and  |xj  indicates  the 
greatest  integer  less  than  or  equal  to  x.  In  Equation  (22),  the  number  of  subcodes  within 
one  code  period  is  Nc  =  nr .  For  example,  if  we  choose  m  =4,  Nc=  16,  the  RNS  folding 
waveform  becomes  Rm=4k=h2 . 16  =  {0,  1,  2,  3,  0,  1,  2,  3,  0,  1,  2,  3,  0,  1,  2,  3}. 


Next,  the  phase  sequences  from  equations  (17)  and  (22)  are  compared.  For 
Equation  (22),  the  parameters  used  to  calculate  the  residue-Frank  phase  code  sequence 
are  m  =  4,  m2  =  16  with  Rm=4Ji=l2_16  =  {0,  1,  2,  3,  0,  1,  2,  3,  0,  1,  2,  3,  0,  1,  2,  3}.  For  the 

Frank  code  from  Equation  (17),  the  parameters  are  M  =  4,  Nc  =  16  with  i,j  =  1,2, 3, 4.  In 
Table  2,  the  phase  code  sequences  are  shown  for  comparison. 

The  phase  code  sequences  of  both  wavefonns  are  identical.  The  ACF,  PACF,  and 
PAF  of  the  CW  waveform  generated  by  Equation  (22)  with  m  =  8,  Nc  =  64  are  also 
identical  to  those  shown  in  Figures  10  and  11,  which  were  generated  by  Equation  (17) 
with  M  =  8,  Nc  =  64. 
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Table  2.  Comparison  of  the  residue-Frank  phase  code  sequence  ( m  =  4  )  with  the 

Frank  phase  code  sequence  (M  =  4,  Nc  =  16  ). 
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D.  RESOLVING  RANGE  AMBIGUITIES  USING  N  RESIDUE-FRANK 
PHASE  CODE  SEQUENCES 

1.  Block  Diagram  of  the  Residue-Frank  Radar  System 


To  extend  the  unambiguous  detection  range,  N  >  2  coprime  moduli  are  used  with 
Equation  (22)  to  generate  the  residue-Frank  phase  code  sequences  necessary  to  modulate 
the  CW  frequency.  The  block  diagram  of  the  radar  using  the  residue-Frank  phase  code 
modulation  is  shown  in  Figure  13.  In  the  transmitter,  the  first  step  is  to  generate  and  store 
the  sequence  values  Rm  k  for  each  modulus.  Then,  a  direct  digital  synthesizer  uses  each 

residue  sequence  Rm  k  to  generate  the  Nc  =  m1  subcodes  (f)m  k  according  to  Equation  (22) 


for  each  modulus.  Each  subcode  corresponding  to  Rmk  has  a  length  of  tb.  Together  the 
Nc  subcodes  have  a  code  period  of  T.  To  detect  the  target’s  unambiguous  range,  each 
phase  code  sequence  must  have  Mms  subcodes.  Each  sequence  is  amplified  and 
transmitted  in  succession. 
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Figure  13.  Block  diagram  for  the  radar  using  N  residue-Frank  phase  code  sequences. 

Upon  reception,  the  signal  is  amplified  and  downconverted  to  an  intennediate 
frequency  compatible  with  an  available  analog-to-digital  (ADC)  converter  technology. 
The  output  of  the  ADC  is  processed  by  a  moving  target  indication  (Elliptic)  filter  to 

remove  the  clutter,  and  then  a  phase  detector  is  used  to  determine  the  phase  and 
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magnitude  of  each  subcode.  The  output  of  the  phase  detector  is  processed  by  a  phase 
code  compressor  which  is  followed  by  a  noncoherent  post-detection  integration  (PDI) 
process.  The  output  of  the  PDI  is  sent  to  a  constant  false  alarm  rate  processor  to  detect 
and  save  the  range  bin  of  the  target  for  that  modulus  sequence.  This  series  of  steps  is 
performed  for  each  residue-Frank  phase  code  sequence,  and  the  detected  target’s  range 
bin  Rm  is  sent  to  a  RNS-to-decimal  algorithm  to  resolve  the  target’s  true  range  after  all 

N  residues  R  are  obtained. 

m 


2.  Transmitted  and  Reference  Codes  for  Compression 

In  this  section,  the  transmission  of  the  set  of  N  residue-Frank  phase  code 
sequences  is  discussed.  The  use  of  each  phase  code  sequence  as  a  reference  sequence  for 
code  compression  is  discussed.  An  example  using  N  =  3  moduli  with  mx  =  3,  m2  =  4, 

and  m3  =  5  is  illustrated.  The  dynamic  range  from  Equation  (21)  for  these  RNS  moduli  is 
=  3x4x5  =  60. 

The  transmitted  sequence  for  each  modulus  requires  the  number  of  subcodes 
equal  to  the  dynamic  range  Mms.  The  transmitted  signal  is  shown  in  Figure  14.  For 

modulus  mt  =3,  there  are  three  residues  within  a  code  period  (0,  1,  2).  For  the  residue- 
Frank  phase  code  period,  the  sequence  has  Nc  =  mf  =  32  =  9.  To  cover  a  dynamic  range 
of  Mms  =  60  subcodes,  the  number  of  code  periods  required  is 


aV  _  ^ RNS 

p,m  2 

m 


(23) 


For  mx  =3,  Np  3  =  60  /  9  =  6.67  code  periods.  The  last  3  subcodes  ( 9  x  7  -  60  =  3  )  are  not 
used  since  they  are  ambiguous. 

For  m2  =  4  and  m3  =  5,  the  number  of  subcodes  Nc  for  each  code  period  are 
ml  =  42  =16  and  m]  =  52  =  25,  respectively.  Also,  the  maximum  number  of  transmitted 
and  received  code  periods  are  Ap4=3.75  (the  last  4  subcodes  are  not  used)  and 
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Np  5  =  2.4  (the  last  15  subcodes  are  not  used),  respectively.  In  Figure  14,  the  code  period 
corresponding  to  each  modulus  is  denoted  as  Tm  where  p  represents  code  period 

number.  Note  that  the  number  of  received  code  periods  processed  in  the  range-Doppler 
matrix  Nr  should  be  less  than  the  number  of  code  periods  received  Np.  That  is 

Nr  <  N  <  N  .  Otherwise,  a  certain  amount  of  correlation  loss  occurs. 

^  Dynamic  range  Mrw  = 60 - H 
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Figure  14.  Illustration  of  transmitted  signal  using  N  =  3  residue-Frank  phase  code 

sequences  for  ml  =  3,  m2  =  4,  and  m3  =  5. 


In  Figure  15  a,  the  residues  Rmk  for  ml  =3,  m2  =4,  and  nu  =5  ( M  RNS  =60)  are 
shown  with  Rmk  =  R  k  2  where  n  e  [0,±1,±2,...| . The  residue-Frank  phase  sequences 
(j>m  k  from  each  modulus  are  shown  in  Figure  15b.  Each  phase  code  sequence  has  a  length 
Nc=m2,  giving  (j)m  k  =  (frm  2 .  If  the  receiver  uses  one  code  period  as  a  reference 
( Nr  =1 ),  the  reference  code  sequence  for  compression  has  a  phase  sequence  from  index 
k  =  1  to  m 2 . 
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(a) 

7"3= 9  m^- 3 


time-index 


(b) 

Figure  15.  Plot  of  the  residue-Frank  phase  code  sequence  for  m1  =  3,  m2  =  4,  and 
m3  =  5  (  M RNS  =  60  )  showing  (a)  residues,  and  (b)  phase  sequences  from  (22). 

In  Figure  16,  the  compression  and  range  bin  detection  process  is  illustrated  for 


mx  =  3  and  Nr  =  1 .  The  residue-Frank  phase  code  sequence  is  returned  to  the  radar 

receiver  with  delay  r  and  then  processed  by  a  filter  matched  to  a  rectangular  pulse  of 
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duration  tb.  The  phase  detector  is  used  to  calculate  the  phase  of  the  signal  which  is 

passed  to  the  compressor.  The  output  from  the  compression  is  then  processed  by  the 
threshold  detection  to  see  whether  the  amplitude  reaches  the  desired  minimum  value.  The 
compressor  processes  a  single  compression  delay  within  a  sampling  period  ts.  If  the 
threshold  detector  detects  the  pulse,  the  presence  of  the  target  will  be  saved  and 
associated  with  residue  Rm.  Each  delay  D,  has  a  duration  of  tb  (e.g.  Dl  is  for  the  time 

delay  between  0  to  tb,  D2  is  for  the  time  delay  between  tb  to  2 tb  and  so  on).  For  m2  =  4 
and  m3  =  5,  the  process  is  similar,  except  the  phase  reference  values  have  to  correspond 
to  the  transmitted  phase  code  sequence  that  is  being  received. 


N— f  — * - T  -  N  th  =  9tb  s*  | 


*  D,  =  fth  range  bin  with  delay  tb  (e.g.  Dl  is  for  r  =  0  to  ,  D2  is  for  r  =  tb  to  2 tb) 


Figure  16.  Illustration  of  compression  at  the  receiver  and  the  range  bin  for  RNS 

ml  =3,  Nc=  9,  and  Nr  =  1 . 
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3.  Calculating  the  Target  Range 

After  the  returned  signal  is  compressed,  the  first  sample  appears  in  the  receiver’s 
range  bin  that  corresponds  to  a  duration  r.  For  each  modular  sequence,  the  target 
detection  corresponds  to  a  residue  Rm.  By  combining  the  paired  residue  from  each 
sequence,  the  target’s  range  bin  can  be  calculated. 

Continuing  the  example  in  the  previous  section,  we  show  the  residues  for 
m1  =  3,  m2  =  4,  and  m3  =  5  and  the  corresponding  target  range  in  Table  3.  Here,  the  signal 

bandwidth  B  =  1  MHz  ,  and  the  subcode  period  th  =  1  /us.  The  unambiguous  range  using 
residue-Frank  phase  code  is 


R 


u,RNS 


RNsh 


(24) 


For  this  example,  the  maximum  unambiguous  range  is  Ru  RNS  =  9  km  (Mms  =  60  )  and  is 

considerably  longer  than  that  of  any  single  code  period.  The  range  resolution  is 
AR  =  ctb  /  2  =  150  m.  The  target’s  range  now  has  to  be  R  <  Ru  ms  in  order  to  be  detected 


with  no  range  ambiguity.  The  range  interval  for  each  range  bin  of  this  example  is  also 
given  in  the  Table  3. 


Table  3.  Bin  matrix  and  range  intervals  of  3  -  channel  RNS  with  ml  =  3,  m2  =  4,  and 

m.  =  5. 
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In  Figure  17,  an  example  of  the  target  detection  process  using  ml  =  3,  m2  =  4,  and 
m,  =  5  is  illustrated.  The  residue-Frank  phase  code  signal  is  transmitted  to  the  target  and 
the  target’s  return  is  in  range  bin  a  =  5  or  770  m.  In  the  code  compression  step,  the 
returned  signal  is  compressed  corresponding  to  each  modulus  in  consecutive  order. 
Shown  in  each  graph  are  residues  Rm  k  (stair-step  waveform),  the  range  bin  index  a  ,  and 
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a  helicopter  target  detection.  The  range  bin  index  a  goes  from  0  to  59  (Mms  =60). 
Since  the  phase  code  sequence  is  periodic,  the  residues  is  also  periodic  with 
R  =  R  ,  2.  For  modulus  m,  =  3,  the  CC  peak  is  repeated  with  T  =  N  t,  =  m:t.  =  9 1,. 

m>k  m,k+nm2  1  5  r  r  mx  c  b  \  b  b 


The  threshold  detector  detects  the  correlated  signal  output,  and  the  range  bin  residue  R, 
is  stored.  Flere,  the  range  bin  residue  corresponding  to  the  correlated  output  for  channels 
m]  =3  is  at  R,  =  2.  The  processes  are  the  same  for  channels  2  and  3  using  the  moduli 

m2  =  4  and  m3  =  5,  respectively.  The  detected  target  range  bin  residues  are  R4  =  1  and 
R5  =  0.  The  residues  are  paired  as  [2  1  0],  and  a  RNS-to-decimal  algorithm  is  used  to 

calculate  the  target’s  true  range  R.  From  Table  3,  the  paired  values  [2  1  0]  correspond  to 
the  range  bin  a  =  5,  and  the  calculated  target  range  is  between  750  m  and  900  m. 


Dynamic 


0 

3 


('fjPiilse  from 
3  ^^2IulJeriod  3 


0  1  2  3  4  5  6  7  8  9  10  11  12  13  14  15  16  17  18  19  20  21  22  23  24  25  26  27  28  29  30 


compression  time  7^  =  vn\ \tb  =  \  6tb 


Qp  ulse  from 
^period 


d 

OI 

•elation 

CO 

ITI 

fi  code 

2 

2o 

ut 

2 

2 

2 

2 

2 

2 

1 

1 

i 

1 

1 

1 

1 

i 

if 

0  1 

T 

I 

£ 

0 

£ 

£ 

_ \JL _ £ 

bin  a 


compression  time  Ts  =  n£tb  =  25 th 


Figure  17.  Example  of  compression  output  for  target  detection  using  N  =  3  residue- 
Frank  phase  code  sequences  with  ml  =  3,  m2  =  4,  and  m3=  5. 
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Although  the  look  up  table  in  Table  3  can  be  used  to  find  the  target’s  range  from 
the  corresponding  residues,  for  a  practical  implementation  the  table  would  be  very  large. 
Alternatively,  the  system  of  simultaneous  congruences  of  the  form  a  =  Rm  (mod  mi )  can 

be  solved  with  a  straightforward  application  of  the  Chinese  remainder  theorem  (CRT). 
The  system  to  be  solved  for  this  example  is  a  =  2(mod3),  a  =  l(mod4),  a  =  0(mod5). 
This  CRT  algorithm  only  requires  the  target  range  bin  residues  Rm  .  The  CRT  solution  of 
the  system  is 


N  M 

7=1  mi 


(25) 


where  a  is  the  range  bin  index.  If  we  consider  this  number  modulo  mi,  we  find  that 


a  =  — ^bRm  =Rm  mod  mr  The  b  values  can  be  found  with  the  Euclidean  algorithm, 

mj  1  1  ’ 

which  is  a  repeated  application  of  the  division  algorithm  [17].  We  note  that  the  7>/  values 

are  also  hardwired  into  the  system  and  depend  only  on  the  moduli  and  not  at  all  on  the 
range  bin  residues. 


4.  Resolving  Multiple  Target  Range  Ambiguities 

In  this  section,  the  detection  of  multiple  targets  is  considered.  As  illustrated  in 
Figure  18,  the  residue  for  the  two  targets  are  periodic  as  Rmk  =  Rm  2 .  If  there  are  no 

missed  targets  in  the  detection  process,  the  sequence  of  residues  within  each  modulus  can 
be  paired  directly  for  the  solution  of  the  target’s  range.  To  solve  for  the  range  of  each 
target,  the  residues  are  paired  as  [2  2  2]  for  the  first  target  and  [2  0  0]  for  the  second 
target.  Using  the  CRT  or  lookup  table,  we  find  the  range  of  the  first  target  is 
150  <  Rl  <  300  m(a  =  2 )  and  the  range  of  the  second  target  is  3000  <  R2  <  3150  m 
(a  =  20). 
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Figure  18.  Resolving  the  range  ambiguity  of  2  targets  using  RNS  for 

m1  =  3,  m2  =  4,  and  m3  =  5. 


In  Figure  19,  if  the  first  set  of  residues  for  the  two  targets  are  not  detected,  the 
residue-Frank  phase  code  waveform  still  allows  the  two  target  range  values  to  be 
determined.  Note  the  residues  are  determined  but  there  is  no  information  available  to  pair 
them  up  correctly.  Additional  information  must  be  used  in  order  to  pair  them  correctly. 
For  example,  for  n?, ,  we  have  R3  =  {2,  2},  for  m2 ,  we  have  R4  =  {2,  0},  and  for  m3,  we 

have  R5  =  {2,  0}.To  pair  them  up  correctly,  the  magnitude  information  associated  with 
each  subcode  (from  the  phase  detector)  can  be  used. 
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Figure  19.  Missed  detections  of  the  first  set  of  residues. 


A  more  convenient  method  is  to  extend  or  fill  in  the  residues  according  to  their 
code  period,  which  can  be  identified  through  examination  of  the  residues  throughout  the 
unambiguous  range.  This  result  gives  the  same  residue  configuration  as  shown  in  Figure 
18.  Note  that  this  method  works  even  with  the  target  residues  overlapping. 


35 


5.  Practical  Considerations 

a.  Detection  of  Target  in  Noise  Using  RNS  Compression 
In  this  section,  the  performance  of  the  N  =  3  residue-Frank  for  mx  =  3, 
m2  =  4,  and  in,  =  5  with  Mms  =  60  is  examined  by  using  MATLAB  to  plot  the  returned 
signals  from  two  targets  at  ranges  of  440  m  and  5840  m.  The  transmitted  signal  has  an 
amplitude  A  =  1,  carrier  frequency  fc  =  1  MFIz,  sampling  frequency  fs=  7  MHz, 

cpp  =  1,  and  tb  =  1  /us.  For  a  practical  case,  the  signals  are  received  with  noise.  The  noise 

changes  the  returned  phase  sequences  and  affects  the  compression  output  value.  Here, 
two  different  SNRs  are  examined.  In  Figure  20,  the  transmitted  signal  without  noise  for 
each  channel  is  plotted.  In  Figure  21,  the  signal  with  SNR  =  30  dB  is  plotted.  The  signal 
power  is  high  compared  to  the  noise  power.  Consequently,  the  noise  signal  does  not 
significantly  affect  the  phase  sequences.  In  Figure  22,  the  SNR  is  reduced  to  0  dB.  The 
phase  sequences  are  changed  significantly. 
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Figure  20.  Residue-Frank  signal  without  noise  for  m,  =  3,  m2  =  4,  and  m3  =  5  with 
Mms  =  60,  A  =  1,  fc  =  1  MHz,  fs  =  7  MHz,  cpp  =  1,  and  th  =  1  ps. 
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Figure  2 1 .  Residue-Frank  signal  with  SNR  =  30  dB  for  ml  =  3,  m2  =  4,  and  m,  =  5 
with  Mms  =  60,  A  =  1,  fc  =  1  MFfz,  fs  =  7  MFIz,  cpp  =  1,  and  tb  =  1  ps. 


Figure  22.  Residue-Frank  signal  with  SNR  =  0  dB  for  ml  =  3,  m2  =  4,  and  m3  =  5  with 
MrnS  =  60,  A  =  1,  fc  =  1  MFIz,  fs=  7  MFIz,  cpp  =  1,  and  tb  =  1  ps. 
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In  Figures  23  and  24,  the  range  bin  profiles  for  each  channel  after 
threshold  detection  of  the  compression  output  for  the  two  targets  are  shown  with  SNR  = 
30  dB  and  SNR  =  0  dB,  respectively.  The  first  target  appears  at  a  range  bin  a  =  2  with 
the  paired  residues  [2  2  2]  (range  300-450  m).  The  second  target  appears  at  range  bin 
a  =  38  with  the  paired  residues  [2  2  3]  (range  5700-5840  m).  Note  that  the  targets 
straddle  the  range  bins  since  the  first  sample  received  from  the  targets  lies  in  the  correct 
range  bin.  Also  shown  is  the  noise  signal  that  is  compressed  before  the  signals  are 
returned  to  the  receiver.  Shown  in  Figure  24,  the  amplitudes  of  compressed  pulses 
degrade  when  the  SNR  is  decreased  to  0  dB.  Also,  the  sidelobes  from  the  compression 
output  are  higher  when  the  noise  power  is  increased.  Another  important  point  is  that  the 
average  sidelobes  from  the  higher  modulus  (e.g.  m3  =  5 )  are  lower  than  those  from  the 

lower  modulus  (e.g.,  ml  =3  and  m2  =4).  This  is  due  to  the  different  processing  gain  or 
Nc  subcodes  that  each  sequence  provides.  For  ml  =3 ,  m2=  4,  and  mi  =  5,  the  processing 
gains  are  mf  =  9,  m2  =  16,  m]  =  25,  respectively. 
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Figure  23.  Target  detection  for  SNR  =  30  dB  showing  the  range  bins  for  m]  =  3,  m2  =  4,  and  m3  =  5  with  the  two  targets  at  440  m 

and  5840  m. 
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Figure  24.  Target  detection  for  SNR  =  0  dB  showing  the  range  bins  for  ml  =  3,  m2  =  4,  and  m3  =  5  with  two  targets  at  440  m  and 

5840  m. 
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One  situation  that  may  occur  is  that  the  threshold  detector  might  read  the 
wrong  range  bin  because  the  amplitudes  of  the  compressed  pulses  fluctuate  due  to  the 
noise  present.  This  results  in  a  target  range  error.  In  Figure  25,  the  first  target  in  channel 
m2  =  4  can  be  seen  with  the  compressed  pulse  detected  by  a  CFAR  (threshold  level  at 
0.8)  at  residue  R4=  3  instead  of  R4  =  2.  So  the  paired  terms  [2  3  2]  (range  7050-7200  m) 

are  computed  instead  of  the  true  range  at  [2  2  2]  (range  300-450  m).  This  is  a  large  range 
error. 


Figure  25.  Illustration  of  detection  error  from  the  threshold  detector  (CFAR)  using 

the  residue -Frank  for  ml  =3 ,  m2=  4,  and  m,  =5. 


b.  Required  CW  Signal  Power  Versus  Required  Output  SNR 

In  this  section,  the  CW  power  required  when  using  the  residue-Frank 
phase  code  sequences  to  detect  the  target  at  the  maximum  unambiguous  range  is 
presented.  Given  the  fixed  unambiguous  range  Ru  =  9000  m,  the  CW  radar  power  can  be 

changed  accordingly  as  a  function  of  the  SNRRo  at  the  output  of  compressor.  The 
relationship  between  the  required  SNRRo  and  the  average  CW  power  ( Pcw  )  is  expressed 
as  [1] 
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(26) 


R:  (4;t)3  kBT0FRBRiLRTLRR  ( SNRRo  t  PG ) 

Cl4/  GtGrA2aTL2 

where  Ru  is  the  unambiguous  range,  kB  =  1 .38 x  ]  0  2  joule/K  (Boltzman’s  constant),  T0 

is  the  ambient  noise  temperature  ( T{)  =  290  K),  FR  is  the  receiver  noise  factor,  BRi  is  the 

radar  receiver’s  input  bandwidth  in  Hz,  Gt  and  Gr  are  the  transmitting  and  receiving 

antenna  gains,  aT  is  the  target’s  radar  cross  section  in  m  ,  L2  is  the  two-way  atmospheric 

transmission  factor,  LRT  is  the  loss  between  the  radar’s  transmitter  and  antenna,  and  LRR 

is  the  loss  between  the  radar’s  antenna  and  receiver. 

In  Figure  26,  the  CW  signal  power  as  a  function  of  the  required  output 
signal-to-noise  ratio  SNRRo  for  ml  =3 ,  m2=  4,  and  m3  =  5  is  plotted  with  BRI  =  1  MHz, 

tb  =  1  jus,  Gt  =  Gr  =  30  dB,  fc=  3  GHz,  A  =  0.1  m,  <jt  =  100  nr,  FR=  5  dB,  L2  =  LRT  = 

Lrk  =  1,  and  Ru  =  cMmsth  12  =  9000  m  (maximum  detection  range  without  ambiguity). 


Figure  26.  Average  power  of  the  CW  transmitter  for  residue-Frank  with 

m1  =  3,  m2  =  4,  and  m3  =5. 
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Note  the  difference  in  the  average  power  required  in  each  channel  due  to 
the  different  PG.  The  channel  m3  =  5,  giving  PG  =  Nc  =  mf  =  25,  requires  less  power 
than  the  other  channels,  ml  =3  (  PG  =  9 )  and  m2  -4  (PG  =  16),  to  detect  the  targets 
within  the  same  unambiguous  detection  range.  For  example,  at  SNRRo  =  30  dB,  the 
required  CW  power  Pcw  for  m3  =  5,  m2  =  4,  and  m,=3  are  7,  10.5,  and  18  W, 
respectively.  Note  that  this  helps  contribute  to  the  LPI  characteristic  of  the  emitter. 

c.  Unambiguous  Detection  Range  for  a  Constant  Output  SNR 

It  is  important  to  see  how  superior  the  residue-Frank  phase  code  is 
compared  to  a  single  Frank  code  with  the  unambiguous  range  Ru  =  cT  /  2  =  cthNc  /  2  and 
limited  to  only  one  code  period.  Here,  the  relationship  between  the  average  signal  power 
Pcw  and  the  maximum  detection  range  Rmnx  is 

R  =  [ - , - - j4 

'  |_(4*)  kJ.FMSNR.J PG)L„L„  J 

However,  the  resolvable  detection  range  is  limited  by  the  unambiguous  range.  The  radar, 
however,  can  detect  the  targets  beyond  Ru  but  there  are  ambiguities. 

In  Figure  27,  the  detection  range  R  =  Rmix  is  plotted  as  a  function  of  the 
required  CW  power  with  constant  SNRRo  =  13  dB,  B  =  1  MHz,  th  =  1  jus,  G,  =  Gr  =  30  dB, 
fc=  3  GHz,  X  =  0.1  m,  aT  =  100  m2,  T0  =  290  K,  FR  =  5  dB,  andC2  =  LRT  =  L M  =1.  Also, 
the  perfonnance  of  each  Frank  code  using  the  corresponding  Nc  is  plotted.  For  the  Frank 
code  (17),  Ml  =  3,  Nc=9  and  M2  =  4,  Nc  =  16, and  M3  =  5,  Nc=25.  Shown  in  the 
graph,  the  single  Frank  code  has  a  limited  unambiguous  range,  where  Ml  =  3  gives 
Ru  =3xl08xl0  6x9/2  =  1350  m,  M2  =  4  gives  Ru  =  3xl08  x  10'6  xl6/ 2  =  2400  m,and 
M3  =  5  gives  Ru  =  3  x  10s  x  10-6  x25  /  2  =  3750  m.  These  values  are  less  than  using  the 
residue-Frank  phase  code,  giving  Ru  =  cMRNSth  /  2  =  9000  m. 
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Figure  27.  Comparison  of  the  maximum  unambiguous  range  of  CW  radar  system  for 
SNRRo  =  1 3  dB  using  the  residue-Frank  phase  code  for  mt  =  3,  m2  =  4,  and 

m3  =  5  with  the  radar  using  each  Frank  code  sequence  individually  with  the 

corresponding  Nc. 

In  conclusion,  the  residue-Frank  phase  code  presented  in  this  chapter  can 
be  used  to  extend  the  unambiguous  range  of  the  CW  polyphase  radar  system.  The  issue 
concerned  with  the  range  bin  error  due  to  the  fluctuation  of  the  compression  outputs  was 
also  introduced.  In  the  next  chapter,  the  relationship  between  the  P4  code  and  the 
symmetrical  number  system  (SNS)  is  presented.  As  will  be  shown  in  examining  the  PAF, 
the  P4  has  a  better  property  in  Doppler  tolerance  than  the  Frank  code.  This  makes  it  more 
attractive  when  the  emitter  attempts  to  detect  the  moving  targets. 
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IV.  P4  POLYPHASE  MODULATION  AND  THE  SYMMETRICAL 

NUMBER  SYSTEM 

In  this  chapter,  the  P4  phase  code  and  the  SNS  are  introduced.  Their  relationship 
is  presented. 

A.  P4  PHASE  CODE 

The  P4  polyphase  code  is  conceptually  derived  from  a  linearly  frequency 
modulated  waveform  (LFMW).  Such  codes  and  compressors  can  be  employed  to  obtain 
much  larger  time-bandwidth  products.  The  significant  advantages  of  P4  codes  are  low 
peak  sidelobes  and  that  they  are  more  Doppler  tolerant  than  other  phase  codes  derived 
from  a  step  approximation  to  an  LFMW  [5],  [12]. 

In  the  P4  code  the  local  oscillator  frequency,  which  is  offset  in  the  /  and  Q 
detectors,  results  in  a  coherent  double  sideband  detection.  The  phase  sequence  of  a  P4 
signal  is  given  by 

(28) 

Ne 

where  k  is  the  time  index  for  a  single  phase  change  per  subcode  and  ranges  from  1  to 
Nc. 

The  phase  values  and  phase  values  modulo  2k  of  a  P4  code  for  Nc  =  64  are 
shown  in  Figures  28  and  29,  respectively.  The  carrier  frequency  is  fc=  1  kHz,  the 
sampling  frequency  is  fs  =  7  kHz,  and  the  number  of  carrier  cycles  per  subcode  is  cpp  = 
1.  The  ACF  and  PACF  for  Nr  =  1  are  shown  in  Figure  30.  The  peak-to-sidelobe  level  is 

PSL  =  20 log10  yj 2  /  ( Nc7i2 )  or  -25  dB.  The  PACF  is  a  perfect  code  with  low  sidelobes. 
The  PAF  in  Figure  31  shows  the  delay-offset  and  Doppler-offset  sidelobes. 
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Figure  28.  P4  phase  sequence  for  Nc  =  64  (From  [1]). 


Figure  29.  Signal  phase  (radians)  modulo  In  versus  k  -index  for  phase  change  of  the 

P4  code  with  Nc  =  64  (From  [1]). 
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(a) 


(b) 


Figure  30.  ACF  (a)  and  PACF  (b)  for  the  P4  code  with  Nc  =  64,  cpp  =  1,  and  Nr  =  1 

(From  [1]). 


Figure  3 1 .  PAF  for  the  P4  code  modulation  with  Nc  =  64,  cpp  =  1,  and  Nr  =  1 

(From  [1]). 

B.  THE  SYMMETRICAL  NUMBER  SYSTEM 

The  SNS  is  composed  of  N  coprime  moduli.  The  integers  within  each  SNS 
modulus  are  derived  from  a  symmetrically  folded  waveform  [18].  The  integer  values 
within  each  SNS  modulus  are  derived  from  a  mid-level  quantization  of  the  symmetrical 
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folding  waveform  and,  therefore,  incongruent  modulus  m  (i.e.,  mod  5:  {0,  1,  2,  2,  1,  0,  1, 
2...}).  Due  to  the  presence  of  ambiguities,  the  set  of  integers  within  each  SNS  modulus 
does  not  form  a  complete  residue  system  by  themselves.  It  is  well  known  that  the 
inclusion  of  additional  redundant  moduli  can  effectively  detect  and  correct  errors  that 
may  occur  within  a  RNS  representation  of  a  number.  The  SNS  formulation  is  based  on  a 
similar  concept,  which  allows  the  ambiguities  that  arise  within  the  SNS  to  be  resolved  by 
using  various  arrangements  of  the  SNS  moduli.  The  vector  of  integer  values  within  a 
single  SNS  folding  waveform  is  given  as 


and 


SL  = 


S„,  = 


0,1,..., 

m 

m 

,-.,2, 1,0,1 

2  J 

2  _ 

.  ,  m  m  ,  _  ,  „  , 

0, 1, —  1, ...,  2, 1, 0, 1 
2  2 


,  when  m  is  odd 


,  when  m  is  even 


(29) 


(30) 


where  [xj  indicates  the  greatest  integer  less  than  or  equal  to  x.  An  example  of  the  SNS 


waveforms  for  ml  =4  and  m2  =  5  (N  =  2)  is  shown  in  Figure  32.  The  period  of  each 
waveform  equals  the  value  of  the  modulus  as 
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Figure  32.  The  SNS  residues  for  ml  =  4  and  m2  =  5. 
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The  dynamic  range  of  the  SNS  M SNS  is 

M sns  =  min  I y I!  mh  +  fl  mi,  \  i22) 

l  1  1=2  l=j+ 1  J 

for  coprime  moduli  mi ,  with  one  of  the  moduli  even,  j  ranges  from  I  to  N  —  l,  N  is  the 
number  of  moduli,  and  mi2,mi3,...,miN  range  over  all  permutations  {2,3,..., N  }.  The 
product  of  j  from  2  to  1  is  empty,  and  its  value  is  1.  If  all  the  coprime  moduli  in  the 
system  are  odd,  then 

f  1  J  l  N  ) 

m  sns = min  t  n  mk + t  n  mi,  (33> 

i^i= i  ^  i= j+ 1  j 

where  j  ranges  from  1  to  N - 1  and  mn,mj2,mi3,...,miN  range  over  all  permutations 
{1,2,3,..., N }.  Note  that  the  dynamic  range  for  the  SNS  is  less  than  the  RNS 

sns  <  ^  RNS  )• 

From  Equation  (32),  for  m1  =  4  and  m2  =  5,  the  dynamic  range  is  M SNS  =  7.  In 
Table  4,  the  two  sequences  are  shown  along  with  the  input  a.  Starting  at  a  =  0,  the  total 
number  of  distinct  paired  terms  is  7,  0  <  a  <  6.  The  index  value  /  =  a  + 1  in  Table  4  is 
shown  to  count  the  integer  number  of  the  distinct  paired  terms  within  the  dynamic  range. 
Also  shaded  is  the  region  with  no  ambiguity. 


Table  4.  Finding  the  SNS  dynamic  range  for  ml  =  4  and  m2  =  5. 
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c. 


SYMMETRICAL  RESIDUE-P4  PHASE  CODE 


The  SNS  symmetrical  residues  Sm  k  can  be  folded  into  the  P4  phase  sequence  by 
associating  Smk  with  the  P4  subcode  phases.  The  phase  index  k  starts  from  1  to  Nc.  The 
symmetrical  residue -P4  phase  sequence  is  given  by 


where  k  =  1,2,3,  ••••,  m.  The  number  of  subcodes  within  one  code  period  is  Nc  =  m.  For 
example,  if  we  choose  m  =  8,Nc=8,  and  the  SNS  sequence  becomes  Sm=sk=l2  g 
=  {0,1, 2, 3, 4, 3, 2,1}. 

Next,  the  phase  sequences  from  equations  (34)  and  (28)  are  compared.  First, 
Equation  (34)  is  used  to  calculate  the  symmetrical  residue -P4  phase  code  sequence  using 
m  =  8,  Nc=  8  with  Sm=s  k=l  ,  8  =  {0, 1, 2, 3, 4, 3, 2, 1} .  For  the  P4  code  in  Equation  (28),  the 

parameters  are  Nc=  8  and  i  =  1,2,3,.  ...8.  In  Table  5,  the  phase  sequences  are  shown  for 
comparison. 

The  phase  code  sequences  of  both  waveforms  are  identical.  The  ACF,  PACF,  and 
PAF  of  the  CW  waveform  generated  by  Equation  (34)  with  m  =  Nc  =  64  are  also 
identical  to  those  shown  in  Figures  30  and  31  which  were  generated  by  Equation  (28) 
with  Nc  =64. 
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Table  5.  Comparison  of  symmetrical  residue-P4  phase  code  sequence 
( m  =  8,  Nc=  8  )  with  the  P4  phase  code  sequence  ( Nc  =  8  ). 
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D.  RESOLVING  RANGE  AMBIGUITIES  USING  N  SYMMETRICAL 

RESIDUE-P4  PHASE  CODE  SEQUENCES 

1.  Block  Diagram  of  the  Symmetrical  Residue-P4  Radar  System 

The  block  diagram  of  a  radar  system  using  symmetrical  residue-P4  phase  code 
sequences  is  shown  in  Figure  33.  This  diagram  is  similar  to  the  residue-Frank  phase  code 
radar  system  in  Chapter  III,  except  that  now  the  symmetrical  residue -P4  phase  code 
sequences  are  used  instead.  Each  phase  sequence  with  the  length  equal  to  the  dynamic 
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range  MSNS  is  transmitted  to  detect  the  target.  In  the  receiver,  each  code  sequence  is 
compressed  with  its  corresponding  conjugate  of  phase  code. 


Symmetrical  Transmit 

residue-Frank  phase  Direct  digital  antenna 


Figure  33.  Block  diagram  for  the  radar  using  N  symmetrical  residue-P4  phase  code 

sequences. 

2.  Transmitted  and  Reference  Codes  for  Compression 

In  this  section,  the  transmission  of  N  symmetrical  residue-P4  phase  code 
sequences  is  explained  as  well  as  the  phase  code  reference  for  compression.  An  example 
using  N  =  3  moduli  with  ml  =  7,  m2  =  8,  and  in,  =  9  is  illustrated.  The  dynamic  range, 

using  Equation  (32)  for  these  SNS  moduli,  is  M SNS  =  37. 

Each  sequence  has  the  number  of  code  length  equal  to  the  dynamic  range  with  the 
transmitted  signal  as  shown  in  Figure  34.  To  create  the  P4  code  sequence,  one  code 
period  requires  the  number  of  subcodes  Nc  =  m.  To  cover  a  dynamic  range  of  subcode, 
the  number  of  code  periods  required  is 

#,,=—■  <35> 

m 

For  mx=l  and  Nc  =  7,  Equation  (35)  gives  N p  7  =  37  /  7  =  5.29,  which  last  5  subcodes 
(7x6-37  =  5)  are  not  included  since  they  are  ambiguous.  For  modulus  m2=  8  and 
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m3  =  9,  the  number  of  subcodes  Nc  for  one  code  period  are  8  and  9  respectively.  Also, 
the  maximum  number  of  transmitted  and  received  code  periods  are  TV  8  =  4.63  (the  last  3 
subcodes  are  not  used)  and  N  9  =  4. 1 1  (the  last  8  subcodes  are  not  used),  respectively.  In 
Figure  34,  the  code  period  corresponding  to  each  modulus  is  denoted  as  Tm  where  p 
represents  the  code  period  number. 


Figure  34.  Illustration  of  transmitted  signal  using  N  =  3  symmetrical  residue -P4 
phase  code  sequences  for  m1  =  7,  m2  =  8,  and  m3  =  9. 


In  Figure  35a,  the  symmetrical  residues  Smk  using  ml  =7,  m2  =8,  and  m3  =  9 
( MSNS  =37  )  are  shown.  Since  the  phase  sequences  are  periodic,  Smk=  Smk+nm  where 
n  e  {0,±1,±2,...} .  The  phase  sequences  from  each  channel  are  plotted  in  Figure  35b.  Each 
symmetrical  residue-P4  phase  sequence  has  a  phase  code  period  equal  to  Nc  =  m,  giving 
<f>m  k  =  (f>m  k+nm .  If  the  receiver  uses  one  code  period  ( Nr  =  1 ),  the  reference  code  sequence 
for  compression  is  the  phase  sequence  from  index  k  =  1  to  in. 
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Figure  35.  Plot  of  the  symmetrical  residue-P4  for  m]  =  7,  m2  =  8,  and  m3  =  9 

(Msns  =  31)  showing  (a)  symmetrical  residues,  and  (b)  phase  sequences  from 

(34). 


In  Figure  36,  the  compression  and  range  bin  detection  process  is  illustrated  for 
m1  =  7  and  Nr  =  1 .  The  symmetrical  residue-P4  phase  code  is  returned  to  the  radar 
receiver  with  delay  z  and  then  processed  by  a  filter  matched  to  the  rectangular  pulse  of 
duration  tb.  Then  the  phase  detector  is  used  to  calculate  the  phase  of  the  signal  and  pass 
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it  through  the  compressor.  The  output  from  the  compression  is  then  processed  by 
threshold  detection  to  see  whether  the  amplitude  reaches  the  desired  minimum  value.  If 
the  threshold  detector  detects  the  pulse,  the  presence  of  the  target  will  be  saved  with  the 
symmetrical  residue  value  Sm.  In  m2  =8  and  =  9,  the  process  is  similar,  except  that 

the  phase  references  have  to  correspond  to  the  transmitted  phase  code  that  is  being 
received. 
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Phase  conjugate 
(reference phase  code) 


Compression  process 
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Figure  36.  Illustration  of  compression  at  the  receiver  and  the  range  bin  for  SNS 

m]  =7,  Nc  =  7,  and  Nr  =  1 . 

3.  Calculating  the  Target  Range 

Using  all  channels,  we  see  that  the  paired  terms  are  mapped  to  a  specific  range  bin 
input  a  (similar  to  residue-Frank  code  bin  detection  process).  Continuing  the  example  in 
the  previous  section  for  ml  =  7,  m2  =  8,  and  m2  =  9,  we  show  the  symmetrical  residue 
paired  values  and  the  corresponding  ranges  in  Table  6.  The  signal  bandwidth  B  =  1  MHz, 
and  the  subcode  period  tb  =  1  /us.  The  unambiguous  range  using  symmetrical  residue-P4 
is 
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(36) 


R 


_  SNS^b 


u,SNS 


For  this  example,  the  maximum  unambiguous  range  is  5550  m  ( MSNS  =37 ).  Range 
resolution  is  1 50  m  ( AR  =  cth  /  2  ). 


Table  6.  Bin  matrix  and  range  intervals  of  3  -  channel  SNS  for  mx  =  7,  m2  =  8,  and 

m3  =  9. 
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In  Figure  37,  an  example  of  target  detection  for  mx  =7,  m2  =8,  and  m,  =9  is 
illustrated.  The  target’s  return  is  in  range  bin  a  =  4  or  650  m.  Shown  in  the  graph  are 
symmetrical  residue  Sm  k  (stair-step  waveform),  the  range  bin  a,  and  a  helicopter  target 

detection.  The  bin  corresponding  to  the  correlated  output  for  channel  m1  =7  is  at 
5  7  =3.  Note  the  repeated  pulse  from  the  compression  output  due  to  the  periodic  phase 
code.  The  process  is  the  same  for  the  channel  m2  =  8,  and  m2  =  9.  The  residue  bins  are 
S  g  =  4  and  S  g  =  4  respectively.  Then  the  three  symmetrical  residues  are  paired  as  [3 

4  4],  and  the  SNS-to-decimal  algorithm  is  used  to  find  the  target’s  true  range  R.  From 
Table  6,  the  paired  terms  [3  4  4]  correspond  to  the  range  bin  a  =  4,  and  the  calculated 
range  is  between  600  m  and  750  m. 


7II2  =  8 


m3=  9 


Figure  37.  Illustration  of  target  detection  by  using  the  range  bin  matrix  for  SNS-P4 

code  with  m[  =  7,  m2  =  8,  and  in,  =  9. 
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In  Chapter  III,  we  solved  the  system  of  congruences  using  CRT  for  recovering  the 
target  range  bin  a  ( <MSNS ).  Due  to  the  definition  of  the  SNS,  a  =  ±Sm  (mod ml ), 
a  =  ±Snh  (mod  m2 ) a  =  ±5  (mod  mN ).  For  each  of  the  congruences,  either  the  plus  or 

minus  is  correct,  but  we  do  not  know  which.  Thus,  we  have  2N  sets  of  N  equations.  For 
example,  if  N  =  2,  we  have 


0 

a  =  Smi  (mod/n, ), 

a  =  Smi  (modm2), 

ii) 

a  =  Smi  (mod  /«, ), 

a  =  -S„h  (mod  m2). 

iii) 

a  =  ~Smi  (mod/n,), 

a  =  ~Sm2  (mod  rn2), 

iv) 

a  =  -Snh  (mod  m, ), 

ci  =  Sm7  (mod/n2  ). 

The  CRT  guarantees  that  each  of  these  has  a  unique  solution  modulo  mim2,  and  exactly 

one  of  these  solutions  lies  within  the  dynamic  range  of  the  system.  This  is  the  value  of  a. 
In  fact,  it  is  only  necessary  to  solve  (i)  and  (ii)  at  most  because  the  solutions  to  (iii)  and 
(iv)  are  the  negatives  of  the  solutions  to  (i)  and  (ii),  respectively. 

Recall  that  in  the  standard  statement  of  the  CRT  we  wish  to  solve  for  a  where 
a  =  Sm  (mod m2),  1  <i<N,  and  the  mi  are  pairwise  relatively  prime.  The  solution  is 


M  DM„b,  a, 


a  —  RNS^\^\  |  ±r*RNSL'2“2  |  _j_  ^_RNS^N^N_ 


m. 


m. 


(37) 


where  Mms  =  ]~j  m; .  Note  that  the  values  of  b;  depend  only  on  mi  and  not  at  all  on 


i=l 


±Sm  .  Thus,  we  may  assume  that  the  constants  ct  =MRNSbi  /  mi  are  known,  and  the  SNS- 
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to-decimal  algorithm  only  needs  to  evaluate  ±q5  ± c2Sm  ±...cNSm^  modulo  Mms  and 
picks  the  one  value  that  lies  within  the  dynamic  range  MSNS. 


4.  Resolving  Multiple  Target  Range  Ambiguities 

In  this  section,  the  detection  of  multiple  targets  is  considered.  In  Figure  38,  we 
see  that  the  symmetrical  residue  for  the  two  targets  are  periodic  as  Sm  k  =  Sm  k+nm.  If  there 

are  no  missed  targets  in  the  detection  process,  the  sequence  of  symmetrical  residues 
within  each  modulus  can  be  paired  directly  for  the  solution  of  the  target’s  range.  To  solve 
for  the  range  of  each  target,  the  residues  are  paired  as  [2  3  4]  for  the  first  target  and  [2  2 
1]  for  the  second  target.  From  the  CRT  or  lookup  table,  the  range  of  the  first  target  is 
750  <  Rx  <  900  m(a  =5  ),  and  the  range  of  the  second  target  is  3900  <  R2  <  4050  m 
(a  =  26). 
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Figure  38.  Resolving  the  range  ambiguity  of  2  targets  using  SNS  for 

m1  =  7,  m2  =  8,  and  m,  =  9. 


In  Figure  39,  if  the  first  set  of  symmetrical  residues  for  the  two  targets  are  not 
detected,  the  symmetrical  residue-P4  phase  code  still  allow  the  two  target  range  values  to 
be  determined.  Note  the  symmetrical  residues  are  determined  but  there  is  no  information 
available  to  pair  them  up  correctly.  Additional  information  must  be  used.  For  example, 
for  mt ,  we  have  S7  =  {2, 2},  for  m2,  we  have  58={3, 2},  and  for  m2,  we  have 
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Sg  =  {4,  1}.  To  pair  them  up  correctly,  the  magnitude  information  associated  with  each 
subcode  (from  the  phase  detector)  can  be  used. 

Rx 

m  i=7 

m  2=8 

m3= 9 


Bin  a 

Figure  39.  Missed  detections  of  the  first  set  of  symmetrical  residues. 

A  more  convenient  method  is  to  extend  or  fill  in  the  symmetrical  residues 
according  to  their  code  periods,  which  can  be  identified  through  examination  of  the 
symmetrical  residues  throughout  the  unambiguous  range.  This  result  gives  the  same 
symmetrical  residue  configuration  as  shown  in  Figure  38.  Note  that  this  method  works 
even  with  the  target  symmetrical  residues  overlapping. 

5.  Practical  Considerations 

a.  Detection  of  Target  in  Noise  Using  SNS  Compression 

In  this  section,  the  performances  of  the  N  =  3  SNS  for  ml  =  7,  m2  =  8,  and 
m3  =  9  with  Msns  =37  are  examined  by  using  MATLAB  to  plot  the  returned  signal  from 
two  targets  at  ranges  of  890  m  and  2990  m,  respectively.  The  transmitted  signal  have 
amplitude  A  =  1, carrier  frequency  fc  =  1  MFIz,  sampling  frequency  fs=l  MHz,  cpp  =  1, 

and  tb  =  1  /us.  The  two  different  SNRs  are  examined.  In  Figure  40,  the  transmitted  signal 
without  noise  for  each  channel  is  plotted.  In  Figure  41,  the  signal  with  SNR  =  30  dB  is 
plotted.  For  the  data  shown  in  Figure  42,  the  SNR  is  reduced  to  0  dB.  The  phase 
sequences  are  changed  significantly  with  the  increase  in  noise  power. 
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Figure  40.  Symmetrical  residue-P4  signal  without  noise  for  m]  =  7,  m2  =  8,  and 
m3  =  9  with  Msns  =  37,  A  =  1,  /'  =  I  MHz,  fs  =  7  MHz,  cpp  =  1,  and  tb  =  1  //s. 
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Figure  41 .  Symmetrical  residue-P4  signal  with  SNR  =  30  dB  for  ml  =  7,  m2  =  8,  and 
m3=  9  with  Msns  =  37,  A  =  1,  fc  =  1  MHz,  fs  =  7  MHz,  cpp  =  1,  and  th  =  1  //s. 
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Figure  42.  Symmetrical  residue-P4  signal  with  SNR  =  0  dB  for  mx  =  7,  m2  =  8,  and 
m3=  9  with  Msns  =  37,  A  =  1,  fc  =  1  MFIz,  fs  =  7  MHz,  cpp  =  1,  and  tb  =  1  //s. 


In  Figures  43  and  44,  the  range  bin  profiles  for  each  channel  are  shown 
using  SNR  =  30  dB  and  SNR  =  0  dB,  respectively.  The  first  target  appears  at  range  bin 
a  =  5  with  the  paired  symmetrical  residues  [2  3  4]  (range  750-900  m).  The  second  target 
appears  at  range  bin  a  =  19  with  the  paired  symmetrical  residues  [2  3  1]  (range  2850- 
3000  m).  Also  shown  is  the  noise  signal  that  is  compressed  before  the  returned  signals 
from  the  targets  arrive  at  the  receiver.  In  Figure  44,  the  amplitudes  of  compressed  pulses 
degrade  when  the  SNR  is  decreased  to  0  dB.  The  effect  from  the  noise  signal  degrades 
the  ability  to  detect  the  targets  similar  to  that  discussed  in  the  residue-Frank  phase  code 
radars.  The  average  sidelobes  from  the  higher  modulus  (e.g.  m3  =9)  are  lower  than  those 

from  the  lower  modulus  (e.g.,  mx  =  7  and  m2=  8)  due  to  the  different  processing  gain. 
For  m1  =  l,m2=  8,  and  my  =  9,  the  processing  gains  are  Ncl  =  7,  Nc2  =  8,  and  Ncl  =  9, 
respectively. 
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Figure  43.  Target  detection  for  SNR  =  30  dB  showing  the  range  bins  for  m]  =  7,  m2  =  8,  and  m,  =  9  with  two  targets  at  890  m  and 

2990  m. 
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Figure  44.  Target  detection  for  SNR  =  0  dB  showing  the  range  bins  for  ml  =  7,  m2  =  8,  and  m3  =  9  with  two  targets  at  890  m  and 

2990  m. 
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Like  the  RNS  phase  code,  a  situation  in  where  the  threshold  detector 
might  read  the  wrong  residue  bin  can  occur  because  the  amplitudes  of  the  compressed 
pulses  fluctuate  due  to  the  noise.  This  results  in  the  target  range  error.  In  Figure  45,  the 
first  target  in  channel  m,  =  8  is  seen  with  the  compressed  pulse  detected  by  CFAR 

(threshold  level  at  0.85)  at  symmetrical  residue  Ss=  2  instead  of  Sg  =3.  So,  the  paired 

terms  are  now  [2  2  4]  (no  range  match)  instead  of  the  true  range  at  [2  3  4]  (range  750- 
900  m).  This  results  in  target  range  error. 


Figure  45.  Illustration  of  detection  error  from  the  threshold  detector  (CFAR)  using 
the  symmetrical  residue-P4  for  mt  =7,  m2=  8,  and  m,  =  9. 

b.  Required  CW  Signal  Power  Versus  Required  Output  SNR 

In  this  section,  the  CW  power  required  when  using  the  symmetrical 
residue-P4  phase  code  sequences  to  detect  the  target  at  the  fixed  unambiguous  range 
Ru  =  5550  m  is  presented.  The  relationship  between  the  required  SNRRo  and  the  average 

CW  power  Pcw  is  given  by  Equation  (26). 

In  Figure  46,  the  CW  signal  power  as  a  function  of  the  required  output 
SNRRo  for  mx  =  7,  m2  =  8,  and  m3  =9  is  plotted  with  BRj  =  1  MHz,  tb  =  1  /us, 
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Gt  =  Gr  =  30  dB,  fc=  3  GHz,  2  =  0. 1  m,  <yT  =  100  m2,  FR=  5  dB,  L2  =  LRT  =  Lm  =  1,  and 
Ru  =  cMSNSth  /  2  =  5550m  (maximum  detection  range  without  ambiguity). 


Figure  46.  Average  power  of  the  CW  transmitter  for  symmetrical  residue-P4  with 

/«,  =  7,  m2  =  8,  and  m3  =  9. 

Note  that  the  difference  of  the  average  power  used  in  each  channel  is  due 
to  the  PG.  The  channel  m3  =  9,  giving  PG  =  Nc  =  mi  =  9,  requires  less  power  than  the 
other  channels,  mx=l  ( PG  =  7  )  and  m2  =  8  ( PG  =  8 ). 

c.  Unambiguous  Detection  Range  for  a  Constant  Output  SNR 

For  a  single  P4  code,  the  unambiguous  range  is  limited  to  only  one  code 
period  (  Ru  =  cT  /  2  =  ctbNc  1 2  ).  In  this  section,  we  compare  the  maximum  unambiguous 
range  using  the  symmetrical  residue-P4  phase  code  sequences  with  the  corresponding 
single  P4  code  sequences.  The  relationship  between  the  average  signal  power  Pcw  and 

the  maximum  detection  range  i?max  is  given  by  Equation  (27). 
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In  Figure  47,  the  detection  range  is  plotted  as  a  function  of  the  required 
CW  power  using  N  =  3  SNS  for  ml  =  7,  m2  =  8,  and  my  =  9  with  SNRRo  =  1 3  dB, 

B  =  1  MHz,  tb  =  1  jus,  G,  =  Gr  =  30  dB,  fc=  3  GHz,  /l  =  0.1  m,  aT  =  100  m2,  T0  =  290  K, 
Fr=  5  dB,and  L2  =  LRT  =  LRR  =1.  For  comparison,  the  performance  of  each  P4  code 

using  the  corresponding  Nc  is  plotted.  For  the  P4  code  from  Equation  (28),  Nci  =  7  and 

Nc2  =  8,  and  Nc2  =  9.  The  single  P4  code  has  a  limited  unambiguous  range,  Ncl  =  7  gives 

Ru  =  3xl08  xl0~6  x7/2  =  1050  m,  Nc2  =  8  gives  Ru  =  3x10s  xl0~6  x8/2  =  1200  m,  and 

Nc3  =  9  gives  Ru  =  3  x10s  xlO-6  x9/2  =  1350  m.  These  values  are  less  than  using  the 

N  =  3  symmetrical  residue-P4  phase  code  sequences,  giving  =  cMSNStb  /  2 

=  5550  m. 


Figure  47.  Comparison  of  the  maximum  unambiguous  range  of  a  CW  radar  system 
for  SNRRo  =  13  dB  using  symmetrical  residue-P4  phase  code  for  ml  =  7, 
m2  =  8,  and  in,  =  9  with  the  radar  using  each  P4  code  individually  with  the 

corresponding  Nc. 
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In  conclusion,  the  symmetric al-P4  phase  code  sequence  presented  in  this  chapter 
can  be  used  to  extend  the  unambiguous  range  of  the  CW  polyphase  radar  systems.  The 
issue  concerned  with  the  range  bin  error  due  to  the  fluctuation  of  the  compression  outputs 
was  considered.  In  the  next  chapter,  the  relationship  between  the  P4  code  and  the  robust 
symmetrical  number  system  (RSNS)  is  presented  to  examine  the  integer  Gray  code 
properties  and  how  they  can  prevent  an  error  in  the  detected  symmetrical  residues. 
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V.  P4  POLYPHASE  MODULATION  AND  THE  ROBUST 
SYMMETRICAL  NUMBER  SYSTEM 

As  discussed  in  Chapters  III  and  IV,  using  the  RNS  and  the  SNS  may  give  a  large 
error  in  the  target’s  range  when  the  threshold  detector  detects  the  wrong  residue  due  to 
the  presence  of  noise  which  can  change  the  phase  of  the  received  signal.  In  this  chapter, 
the  RSNS  is  presented,  and  it  is  shown  that  the  range  errors  in  the  robust  symmetrical 
residue-P4  phase  code  sequence  can  be  bounded  to  a  small  value  due  to  the  Gray  code 
property  of  the  RSNS  [20].  The  P4  phase  expression  is  presented  and  examples  are  given 
to  illustrate  the  target  detection  process  with  this  new  waveform. 

A.  P4  PHASE  CODE 

The  details  of  the  P4  phase  code  are  discussed  in  Chapter  IV.  The  P4  phase  code 
is  derived  from  a  linearly  frequency  modulated  waveform  (LFMW).  To  calculate  the 
phase  sequence,  Equation  (28)  is  used.  Recall  the  characteristics  of  the  P4  code  sequence 
is  its  superior  PAF  sidelobe  structure  and  its  Doppler  tolerance  properties. 

B.  THE  ROBUST  SYMMETRICAL  NUMBER  SYSTEM 

The  RSNS  is  a  modular  system  consisting  of  N  >  2  integer  sequences  with  each 
sequence  associated  with  a  coprime  modulus  mi .  Due  to  the  presence  of  ambiguities,  the 
set  of  integers  within  each  RSNS  sequence  do  not  form  a  complete  residue  system.  A  set 
of  integers  y0 ,  y, , . . . ,  ym_x  (mod  m)  form  a  complete  residue  system  if  y  =  g(modw)  for 

g  =  0,l,...,m-l  [21].  The  ambiguities  within  each  modulus  sequence  are  resolved  by  the 
use  of  additional  moduli  and  considering  the  vector  of  paired  integers  from  all  N 
sequences.  The  RSNS  is  based  on  the  sequence 

[RS'm  h }  =  [0, 1, 2, . . . ,  m  - 1,  m, m  - 1, ....  2, 1]  (38) 

for  an  integer  h  such  that  0  <  h  <  2m.  To  form  the  A-sequence  RSNS,  each  term  in  (38) 
is  repeated  N  times  in  succession.  The  integers  within  one  folding  period  of  a  sequence 
are  then 
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(39) 


m,m,  •  ■  ■  m,  m  - 1, •  •  • ,  m  - 1,  •  •  • ,  1,  •  •  • ,  1]. 

This  results  in  a  periodic  sequence  with  a  period  of  [20] 

Pm  =  2Nm.  (40) 

That  is,  each  RSNS  sequence  contains  mi  integers  with  the  folding  period  of  each 
sequence  equal  to  2 Nmi.  The  construction  of  the  N  sequences  ensures  that  any  two 

successive  RSNS  vectors  (paired  integers  from  the  N  sequences)  differ  by  only  one 
integer,  resulting  in  an  integer  Gray  code  property. 

Each  sequence  corresponding  to  mi  is  right  (or  left)  shifted  by  st  =  i  - 1  places. 
The  chosen  shift  values  {  Sj ,  s2 , . . . ,  sN  }  must  also  form  a  complete  residue  system  modulo 
N.  Further,  if  each  sequence  is  extended  periodically  with  period  2 Nm  as 
RSmM„2Nm  =  RR,„j,  where  n  e  |0,±1,±2,...},  then  RSmh  is  called  a  symmetrical  residue  of 
( h  +  n2Nm  )  modulo  2 Nm. 

The  calculation  of  M RSNS  and  its  location  in  the  sequence  are  a  function  of  N  and 

the  choices  for  mi  and  st .  A  closed-form  solution  for  computing  MRSNS  for  N  -2 
moduli  is  reported  in  [22],  and  summarized  as  follows: 


For  ml  >  3  and  m2  =  ml  + 1 


M rsns  =  3(/h,  +  m2 )  -  6  =  6  ml  -  3 . 

(41) 

For  5  <  ml  <  m2  and  m2  <  m]  +  2 

M RSNS  ~  ^  • 

(42) 

For  5  <  m1  <  m:  and  m2  >  mt  +  3 

M rsns  ~  2in2  —  2. 

A  closed-form  solution  for  TV  =  3  moduli  of  the  form  2~  -1,2’,  2"  +1  is 

y  3  2  15 

M rsns  —  2  ~l~  2  m'  +  ^ 


(43) 


(44) 
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where  z  is  any  integer  and  ml  >  3  [23],  An  efficient  algorithm  for  computing  M RSNS  and 
its  position  for  any  general  set  of  moduli  is  reported  in  [24]— [25] . 

The  shift  values  s.  for  each  sequence  do  not  affect  the  RSNS  dynamic  range  but 
do  make  a  difference  in  the  location  of  the  beginning  position  h  of  the  dynamic  range 
and  the  ending  position  h  +  MRSNS  - 1 .  After  choosing  the  beginning  h  ,  it  is  useful  to 
know  the  symmetrical  residues  at  this  point  (e.g.,  to  align  the  symmetrical  residues  for 
each  sequence).  Since  each  folding  period  consists  of  2 Nmi  integers,  the  symmetrical 
residues  are  determined  by  first  subtracting  off  an  integer  number  of  2 Nmi  integers  as 


n,  =  h  - 


2  Nm. 


2  Nm.  . 


(45) 


This  value  is  then  used  to  find  the  symmetrical  residue  RSm  h  as  [20] 


RSm,h=\ 


N 


2 Nmi  +  N -ni  +  si  - 1 
N 


st  <  ni  <  Nmi  +  si  + 1 
Nm..  +  s..  +  2  <  n,  <  2 Nm.  +  v  - 1 . 


(46) 


Let  Xh  be  the  vector  of  N  paired  integers  from  each  sequence  at  h  where  h  is  the 
position  of  the  vector  X h  within  the  RSNS.  In  Table  7,  h  and  X h  for  an  N  =  3  RSNS 
system  with  mj  =[3  4  5]r  and  right  shift  sf  =[0,l,2]r  are  shown.  For  this  example,  the 
RSNS  dynamic  range  M RSNS  =43,  and  the  position  begins  at  /?  =  61  [20].  The  vector 

*61  =  [2  4  I  ]7 .  The  set  of  integers  that  lie  within  the  dynamic  range  MRSNS  =43  contain 
no  ambiguities.  Also  note  in  Table  7  the  integer  Gray  code  property  where  any  code 
transition  results  in  just  one  integer  changing  value  by  ±1 .  Also  shown  is  the  bin  index  a 
which  runs  from  0  to  42  (total  43)  that  will  be  used  to  refer  to  the  range  bin  detection. 
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Table  7.  The  RSNS  folding  waveforms  for  m]  =  3  (sl  =  0),  m2  =  4  (s2  =  1),  and 

m,  =5  (s3  =  2)  (After  [20]). 
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C.  ROBUST  SYMMETRICAL  RESIDUE-P4  PHASE  CODE 

Similar  to  the  SNS  phase  code,  the  RSNS  folding  wavefonn  can  be  related  to  the 
P4  code  sequence.  The  RSNS  symmetrical  residues  RSmk,  can  be  folded  into  the  P4 

phase  sequence  by  associating  RSm  k,  with  the  P4  subcode  phases.  The  phase  index  k' 
starts  from  1  to  A  ,  where  A  is  the  number  of  RSNS  subcodes  within  one  code 
period.  The  robust  symmetrical  residue-P4  phase  sequence  is  given  by 


fim.k'  ~ 


n 
2  m 


(RSm,k'  ~m)2 


n 

- m 

2 


and 


A,  =  2  mN 

CRSNS 


(47) 


(48) 


where  A  is  the  number  of  channels.  For  example,  if  we  choose  m  =  4,  A  =  2,  and 
N  =  16.  The  robust  symmetrical  residue-P4  phase  code  sequence  is  RSm=4k,= 
{0,0, 1,1, 2, 2, 3, 3, 4, 4, 3, 3, 2, 2, 1,1}.  Note  that  the  folding  sequence  is  similar  to  SNS  with 
m  =  8,  Sm=sk  =  {0,1, 2, 3, 4, 3, 2,1}  except  that  the  robust  symmetrical  residue-P4  phase  is 

repeated  A  times;  there  are  A  phase  values  for  each  robust  symmetrical-P4  phase  for 
every  P4  phase  value.  That  is,  both  k'  =  1  and  k'  =  2  of  the  RSNS  corresponds  to  k  =  1 
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for  the  P4  phase.  Consequently,  the  number  of  subcodes  Nc  for  the  P4  code  is  related  to 
the  number  of  RSNS  subcodes  N  as 

CRSNS 

Nc 

N  =  — (49) 
e  N 

Also,  the  PG  for  the  robust  symmetrical  residue-P4  waveform  is 

PGRSNS=Nc=^  =  2m.  (50) 

1.  Example  of  Phase  Sequence  Calculation 

In  this  section,  the  phase  code  sequences  from  equations  (47)  and  (28)  are 
compared.  Equation  (47)  is  used  to  calculate  the  phase  sequence  for  m  =  4,  N  =  2,  and 

N  =  16.  Using  the  RSNS,  we  find  that  the  robust  symmetrical  residues  are 

RSm=Ak.  =  {0,  0,  1,  1,  2,  2,  3,  3,  4,  4,  3,  3,  2,  2,  1,  1}  where  k'  =  1,2,..., 16.  In  order  to 

compare  to  the  P4  phase  sequence  in  Equation  (28)  with  Nc  =  NCrns  / N  =  16/2  =  8,  k  = 

1,2, 3,... ,8,  the  two  phase  sequences  are  shown  in  Table  8. 


Table  8.  Comparison  of  the  robust  symmetrical  residue-P4  phase  code  sequence 
( m  =  4,  N  =  2,  Nr  =16)  with  the  P4  phase  code  (  N=  8). 

X  CRSm  /  A  V  C  ' 


Phase 

index 

(k') 

Robust  symmetrical  reside-P4  phase  code 

(47) 

P4  phase  code  (28) 

aw 

<t>m= 4,A-'(rad) 

k 

<t>k  (rad) 

1 

0 

1 

O 

2 

0 

^=CiU0-4>2)-f4-0 

3 

1 

^=[S(1“4)2]“f4=_2'7489 

2 

2 

^2  =  ^2-i)  n{2  1)=  2.7489 

8 

4 

1 

^4=(S(1“4)2]“f4=_2'7489 
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5 

2 

<h=[ 

^2-4'2) 

4=-4.7124 

2 

3 

2 

=  ,r(3  1)=  4.7124 

8 

6 

2 

A(2-4>2) 

4=-4.7124 

2 

7 

3 

<h=[ 

A3-4)2) 

4=-5.8905 

2 

4 

^4  =  ;r(4_1)  ;r(4  1)=  5.8905 

8 

8 

3 

<k=\ 

>4>2) 

4=-5.8905 

2 

9 

4 

<h=[ 

;^4-4>2) 

4=-6.2832 

2 

5 

2 

^5  =  /r(5-D  ^  1)=  62g32 

8 

10 

4 

<ho= 

|-C4=_6.2832 

11 

3 

<h\= 

I-— 4=-5.8905 

1  2 

6 

2 

^6=^-(6-d  ^(6  1)=  5  8905 

8 

12 

3 

<h2= 

4=-5.8905 

1  2 

13 

2 

f^3  = 

(s<2-4>2: 

— 4= — 4.7 1 24 

1  2 

7 

<h  =  K(7~lY  ^(l-l)=-4.7UA 

O 

14 

2 

j-y4=-4.7124 

15 

1 

<k5= 

I-— 4=-2.7489 

1  2 

8 

2 

^  =  ^(8-D-_;r(8_i)=_2.7489 

O 

16 

1 

<h6= 

[U'-'A 

4=-2.7489 

2 

2.  ACF,  PACF,  and  PAF 

In  this  section,  we  examine  the  ACF,  PACF,  and  PAF  of  the  signal.  To  compare 
the  robust  symmetrical  reside-P4  phase  code  sequence  with  the  P4  phase  code  sequence, 
the  ACF,  PACF,  and  PAF  are  examined  for  each  waveform.  In  Figure  48,  the  robust 
symmetrical  residues  using  mi  =  16,  A  =  2,  and  NCg  ^  =64  are  shown.  In  Figure  49,  the 

robust  symmetrical  residue -P4  phase  code  sequence  residues  are  shown  (using  Equation 
(47))  with  fc  =  1  kHz,  fs  =  7  kHz,  and  cpp  =  2.  Note  that  cpp  represents  the  number  of 
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carrier  cycles  per  phase  change,  and  the  cpp  value  for  the  P4  is  N  times  that  of  the 
robust  symmetrical  residue-P4  phase  code.  For  the  P4,  the  subcode  period  is 


h  ~ 


cpp 

fc  ' 


For  the  robust  symmetrical  residue-P4,  the  subcode  period  is 


(51) 


(52) 


From  (5 1)  and  (52),  we  have 

th=Nt'h.  (53) 

Using  Equation  (48),  we  see  that  the  phase  sequence  from  the  robust  symmetrical 
residue-P4  for  m  =16,  N  -  2,  and  N  =  64  is  actually  the  same  as  the  P4  phase  code 

1  CRSNS  J  r 


sequence  using  Nc  =32. 


Figure  48.  RSNS  residues  RSm  ,  for  modulus  m  =  16,  N  =  2,  and  Nr  =  64. 

b  m>k  ?  5  Crsns 
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Figure  49.  The  robust  symmetrical  residue-P4  phase  code  sequence  using  m  =  16, 

N  =  2,  NCr^s  =  64  is  shown  with  the  relationship  with  the  P4  code  using  Nc  =  32. 

The  power  spectral  density  (PSD)  of  the  signal  is  shown  in  Figure  50.  The  signal 
bandwidth  B  is  reduced  by  the  factor  of  N.  With  cpp  =  2,  the  bandwidth  of  the  robust 
symmetrical-P4  signal  is  B  =  fcl  cpp  =  1  /  tb  =  500  Hz  (instead  of  B  =  1  kHz  when 
cpp  =  1).  In  Figure  51,  the  ACF  and  PACF  are  plotted.  The  signal’s  characteristics  are 
the  same  as  the  P4  phase  code  sequence  using  Nc  =  32,  th  =  2  ms.  In  the  ACF  plot,  the 

output  gives  PSL  =  201og10 1^2/  (Nc7T2)  j  =  -22  dB.  Note  that  the  width  of  the  mainlobe 

in  the  ACF  plot  is  2 t'b  and  is  proportional  to  the  number  of  moduli  N.  That  is,  the  higher 
the  number  of  moduli  A,  the  broader  the  mainlobe.  In  Figure  52,  the  PAF  is  plotted.  The 
PAF  from  the  robust  symmetrical  residues  using  mi=\6,N  =  2,  and  Ar^  =  64  in 
Figure  52a  is  compared  to  the  PAF  from  the  P4  code  sequence  using  Nc  =  32  in  Figure 
52b.  Both  PAFs  are  the  same. 
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Figure  50. 
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Power  spectral  density  of  the  robust  symmetrical  residue-P4  signal  using 
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Figure  5 1 .  The  ACF  and  PACF  of  RSNS  phase-coded  signal 

( m  =  16,  N  =  2,  N  =  64  )  at  the  receiver  with  N  =  1, 

v  CRSNS  y  r 

fc  =  1  kHz,  fs  =  7  kHz,  and  cpp  =  2. 
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^  NcRSNStb 


(a) 


(b) 

Figure  52.  The  PAF  of  the  signal  for  (a)  the  robust  symmetrical  residue-P4  phase 
code  sequence  using  m  =  16,  N  =  2,  N  =  64,  and  t'h  =  1  ms,  and  (b)  the  P4  code 

using  Nc  =  32,  and  tb=  2  ms. 
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D.  RESOLVING  RANGE  AMBIGUITIES  USING  N  ROBUST 
SYMMETRICAL  RESIDUE-P4  PHASE  CODE  SEQUENCES 

1.  Block  Diagram  of  the  Robust  Symmetrical  Residue-P4  Radar  System 

To  extend  the  unambiguous  detection  range,  N  >  2  coprime  RSNS  sequences  are 
used.  The  block  diagram  of  the  radar  using  A -channel  RSNS  phase  coding  is  shown  in 
Figure  53.  This  diagram  is  similar  to  the  residue-Frank  phase  code  radar  system  except 
that  now  the  RSNS  moduli  are  used  instead.  Each  phase  sequence  with  the  length  equal 

to  the  dynamic  range  MRSNS  is  transmitted  to  detect  the  target.  In  the  receiver,  each  code 
sequence  is  compressed  with  its  corresponding  reference  phase  code. 
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Figure  53.  Block  diagram  for  the  radar  using  N  robust  symmetrical  residue-P4  phase 

code  sequences. 


2.  Transmitted  and  Reference  Codes  for  Compression 

In  this  section,  the  transmission  of  a  set  of  N  robust  symmetrical  residue-P4 
phase  coded  sequences  is  explained.  Also,  the  use  of  the  phase  code  sequences  as  the 
reference  sequence  for  code  compression  is  discussed.  An  example  using  N  =  3  RSNS 
moduli  with  ml  =  3,  m2  -  4,  and  m,=  5  is  illustrated.  The  dynamic  range  for  these  RSNS 

moduli  from  Equation  (44)  is  MRSNS  =43.  The  starting  position  h  and  ending  position 
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h+M RSNS  -1  are  61  and  103,  respectively  [20],  The  robust  symmetrical  residues  within 
the  RSNS  moduli  are  paired  as  shown  in  Table  7. 

The  transmitted  signal  is  shown  in  Figure  54.  Each  sequence  has  a  length  equal  to 
the  dynamic  range.  For  ml  =  3,  there  are  four  robust  symmetrical  residues  within  the 

sequence,  RS3  =0,  1,  2,  and  3.  To  create  the  P4  code  sequence,  one  code  period  needs 
the  number  of  subscodes  N  =2 Nm,  =  2x3x3  =  18.  For  m,  =  4,  the  robust 

symmetrical  residues  are  0,  1,  2,  3,  and  4.  The  required  number  of  subcodes  is 
2x3x4  =  24.  For  m3  =5,  the  robust  symmetrical  residues  are  0,  1,  2,  3,  4,  and  5.  The 
number  of  subcodes  is  2x3x5  =30.  As  shown  in  the  Figure  54,  the  beginning  of  each 
sequence  is  not  necessarily  the  beginning  of  the  phase  code  sequence  k’  =  1  (e.g.,  the 
beginning  phase  sequences  for  m,  =3,  m2  =  4,  and  m3  =  5  are  k'  =  8,  k'  =  13,  and 
k'  =  30,  respectively.  The  reference  code  sequences  in  the  compressors  can  be  chosen  to 
be  either  starting  at  the  beginning  of  the  sending  signal  waveforms  or  the  starting  of 
phase  sequence  at  k'  =  1 .  The  importance  is  that  the  reference  code  has  to  have  a  length 
equal  to  the  code  period  Tm  for  each  channel  in  order  to  correctly  calculate  the  code 
compression  output.  In  this  example,  we  choose  to  use  the  reference  signals  starting  from 
the  beginning  of  the  transmitted  signal  ( 0S  to  07  for  m]  =  3,  013  to  012  for  m2  =  4,  and 

03o  to  0,9  for  m3  =  5  ).  The  number  of  sequence  periods  is  calculated  as 

A  _  M RSNS  _  Mrsns  rszn 

p’m  Nc  2 mN  '  1  ' 

CRSNS 

Sending  43  subcodes,  the  modulus  m]  =  3  needs  N ,  =43/18  =  2.4 .  For  modulus 
m2  =  4  and  m3  =  5,  the  number  of  subcodes  N  are  24  and  30,  respectively.  Also,  the 
number  of  transmitted  code  periods  are  Np4  =1.8  and  Np5  =1.4,  respectively.  Here, 
only  channel  m1  =  3  is  sending  more  than  one  complete  code  period.  So  we  can  expect  to 
see  two  compressed  pulses  from  channel  m1  =3  if  all  subcodes  return  back  to  the 
receiver  (there  will  be  only  one  pulse  in  the  other  two  channels). 
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Figure  54. 


Illustration  of  transmitted  signal  using  A  =  3  robust  symmetrical  residue- 
P4  phase  code  sequences  for  m]  =  3,  m2  =  4,  and  m,  =  5. 


In  Figure  55a,  the  RSNS  folding  waveforms  for  777,  =  3,  m2  =  4,  and  m3  =  5  with 
M R SNS  =  43  are  plotted.  Also,  the  phase  sequences  from  each  channel  are  plotted  in 
Figure  55b.  Since  RS mMn2Nm  =  RS mh  where  n  e  {0,±l,±2,...},the  phase  sequences  are 
periodic  as  </>mX+n2Nm  =  <j)mX . 

In  Figure  56,  the  compression  and  range  bin  detection  process  is  illustrated  for 
RSNS  m1  =3,  and  Nr  =  1.  The  returned  phase  coded  signal  is  first  processed  by  a  filter 

matched  to  the  rectangular  pulse  of  duration  t'b.  The  phase  detector  is  used  to  calculate 
the  phase  of  each  sample  and  then  pass  it  through  the  compressor.  The  output  from  the 
compression  is  then  processed  by  the  threshold  detection  to  see  whether  the  amplitude 
reaches  the  desired  minimum  value.  If  the  threshold  detector  detects  the  pulse,  the 
presence  of  the  target  will  be  saved  with  the  robust  symmetrical  residue  value  RSm.ln 

m2=  4  and  m3  =  5,  the  process  is  similar,  except  that  the  phase  references  have  to  match 
with  the  transmitted  phase  code  that  is  being  received. 
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Figure  55. 
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Figure  56.  Illustration  of  compression  at  the  receiver  and  the  range  bin  for  RSNS 

m1  =3,  Nc  =  18,  and  Nr  =  1 . 


3.  Calculating  the  Target  Range 


Using  all  sequences,  we  map  the  paired  terms  to  a  specific  range  bin  a  (similar  to 
the  residue-Frank  code  bin  process  in  Chapter  III).  Continuing  the  example  in  the 
previous  section  for  mx  =3 ,  m2=  4,  and  m3  =  5,  we  show  the  robust  symmetrical  residue 
paired  values  and  the  corresponding  ranges  in  Table  9.  The  signal  bandwidth  is 
B  =  1  MHz,  and  the  subcode  period  is  t'h=tb  /  N  =  0.33  /us.  The  new  expression  formula 
for  calculating  the  unambiguous  range  using  RSNS  is 


R 


u, RSNS 


cA4  t' 

RSNSlb 


(55) 


For  this  example,  the  maximum  unambiguous  range  is  Ru  RSNS  =  2150  m.  The  range 


resolution  is  50  m  (  AR  =  ct'b  /  2  ). 
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Table  9.  Bin  matrix  and  range  intervals  of  3  -  channel  RSNS  for  ml  =  3,  m2  =  4,  and 

m3  =  5. 
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In  Figure  57,  an  example  of  target  detection  for  ml  =3,  m2=  4,  and  m3  =  5  is 
illustrated.  The  return  signal  is  in  the  range  bin  a  =  5  or  270  m.  First,  the  phase  code 
signal  of  the  modulus  ml  =  3  (top)  is  compressed.  Recall  that  the  RSNS  has  a  width  of 
the  compression  output  that  covers  N  =  3  code  sequences  at  robust  symmetrical  residues 
RS,  =  3,  2,  and  2.  So  the  threshold  detection  has  to  decide  to  take  one  bin,  e.g.,  at  the 
peak  of  the  compressed  pulse,  giving  the  robust  symmetrical  residue  RS3  =  2.  The 
process  is  the  same  for  the  channels  m2  =  4  and  m3  =  5.  The  robust  symmetrical  residues 
are  RS4  =  3  and  RSs  =  1,  respectively.  The  three  robust  symmetrical  residues  are  paired 
as  [2  3  1]  and  must  be  converted  to  the  decimal  RSNS  bin  number.  In  Table  9,  the 
corresponding  range  is  between 250  m  and  300  m  (bin  a  =  5). 


Figure  57.  Illustration  of  target  detection  by  using  the  range  bin  matrix  for  RSNS-P4 

code  with  ml  =3,  m2  =  4,  and  m3  =  5. 
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4.  Resolving  Multiple  Target  Range  Ambiguities 

In  Figure  58  for  m]  =3,  the  compressed  pulses  are  repeated  as  the  returned  code 
runs  through  the  reference  code  in  the  compressor.  The  two  target  returns  overlap  one 
another  at  range  bins  a  =  19,  20,  and  21  and  are  ambiguous.  Next  we  resolve  the  range 
ambiguity  by  using  3  channels.  If  the  threshold  detector  is  set  to  measure  at  the  peak  of 
the  compression  output,  the  robust  symmetrical  residues  from  the  first  target  are  paired  as 
[3  4  0],  giving  a  range  bin  a  =  2  (range  50-100  m).  Similarly,  the  robust  symmetrical 
residue  paired  terms  from  the  second  target  are  [3  2  4],  giving  a  range  bin  a  =  20  (range 
1,000-1,050  m).  The  two  targets  have  different  paired  terms,  so  the  range  ambiguity  can 
be  resolved. 
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Figure  58.  Resolving  the  range  ambiguity  of  2  targets  using  RSNS  for 

m1  =  3,  m2  =  4,  and  m3  =  5. 


As  mentioned  in  Chapters  III  and  IV,  we  can  resolve  the  targets  even  if  the  first 
set  of  the  residues  are  not  detected.  Imagine  we  use  a  higher  set  of  RSNS  moduli  and  the 
dynamic  range  is  higher  with  the  result  that  the  code  period  repeats.  To  find  the  true 
range  of  the  targets  that  are  missed,  the  magnitude  information  or  the  filling  in  the  range 
bin  according  to  the  code  periodicity  can  be  used. 
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5.  Practical  Considerations 

a.  Detection  of  Target  in  Noise  Using  RSNS  Compression 

In  this  section,  the  performances  of  the  N  =  3  RSNS  for  ml  =  3,  m2  =  4, 

and  m3  =  5  with  M  RSNS  =  43  are  examined  by  using  MATLAB  to  plot  the  returned  signal 
from  two  targets  at  ranges  of  140  m  and  1040  m.  The  transmitted  signal  has  amplitude 
A  =  1,  carrier  frequency  fc  =  1  MHz,  sampling  frequency  fs=  9  MHz,  cpp  =  1,  th  =  1  /us, 

and  t'h  =  0.33  /us.  In  Figure  59,  the  transmitted  signal  without  noise  for  each  channel  is 

plotted.  In  Figure  60,  the  signal  with  SNR  =  30  dB  is  plotted.  The  signal  power  is  high 
compare  to  the  noise  power.  As  a  result,  the  noise  signal  does  not  significantly  affect  the 
phase  sequences.  In  Figure  61,  the  SNR  is  reduced  to  0  dB.  The  phase  sequences  are 
changed  significantly. 
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59. 

Robust  symmetrical  residue-P4  signal  without  noise  for 

,  m2 

=  4,  and  m,  =  5  with  M  RSNS  =43,  A-  1,  fc 

=  1  MHz,  f  =9  MHz, 

cpp 

=  1,  and  t'b  =  0.33  /us. 
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Figure  60.  Robust  symmetrical  residue-P4  signal  with  SNR  =  30  dB  for 
ml  =  3,  m2  =  4,  and  m,  =  5  with  M RSNS  =43,  A  =  1,  fc  =  1  MHz,  fs  =  9  MHz, 

cpp  =  1,  and  t'b  =  0.33  ps. 
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Figure  61.  Robust  symmetrical  residue -P4  signal  with  SNR  =  0  dB  for 
m,  =  3,  m2  =  4,  and  m3  =  5  with  MRSNS  =43,  A  =  1,  fc  =  1  MHz,  fs  =  9  MHz, 

cpp  =  1,  and  t'b  =  0.33  ps. 
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In  Figures  62  and  63,  the  range  bin  profiles  from  two  targets  are  shown 
using  SNR  =  30  dB  and  SNR  =  0  dB,  respectively.  The  first  target  appears  at  range  bin 
a  -  2  with  the  paired  robust  symmetrical  residues  [3  4  0]  (range  100-150  m).  The  second 
target  appears  at  range  bin  a  =  20,  giving  paired  terms  [3  2  4]  (range  1000-1050  m).  As 
shown  in  Figure  63,  the  amplitudes  of  the  compressed  pulses  degrade  when  the  SNR  is 
decreased  to  0  dB.  Also,  the  sidelobes  from  the  compression  output  are  higher  when  the 
noise  power  is  increased.  Another  important  point  is  that  the  average  sidelobes  from  the 
higher  modulus  (e.g.  m3  =5)  are  lower  than  those  from  the  lower  modulus  (e.g.  m]  =  3 

and  m2  =4)  due  to  the  different  processing  gains.  For  m]  =3,  m2  =4,  and  in,  =5,  the 
processing  gains  are  Ncl  =  2mx  =  6,  Nc2  =  2m2  =  8,  and  AN  =  2m,  =  1 0.  respectively. 
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Figure  62.  Target  detection  for  SNR  =  30  dB  showing  the  range  bins  for  ml  =  3,  m2  =  4,  and  m3  =  5  with  two  targets  at  140  m  and 

1040  m. 
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Figure  63.  Target  detection  for  SNR  =  0  dB  showing  the  range  bins  for  m]  =  3,  m2  =  4,  and  m3  =  5  with  two  targets  at  140  m  and 

1040  m. 
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Using  N  robust  symmetrical  residue-P4  phase  code  sequences  allows 
control  of  the  range  coding  error.  The  threshold  detector  might  read  the  wrong  robust 
symmetrical  residue  bin  because  the  amplitudes  of  the  compressed  pulses  are  fluctuating. 
However,  the  paired  terms  from  the  robust  symmetrical  residues  will  change  only  one 
position  (the  Gray  code  property).  In  Figure  64,  the  compressed  pulse  from  the  first  target 
in  m2  =  4  is  detected  by  CFAR  (threshold  level  at  0.85)  at  the  robust  symmetrical  residue 

RS4  =  3  instead  of  RS4  =  4.  The  paired  terms  are  now  [3  3  0]  (range  150-200  m)  instead 

of  the  true  range  at  [3  4  0]  (range  100-150  m).  The  maximum  range  error  is  only  one 
range  bin  or  50  m. 


Figure  64.  Illustration  of  detection  error  from  the  threshold  detector  (CFAR)  using 
the  robust  symmetrical  residue-P4  for  m]  =  3,  m2  =  4,  and  m3  =  5. 


b.  Plot  of  Required  Signal  Power  Versus  Required  SNR  Output 

In  this  section,  the  power  required  when  using  the  robust  symmetrical 
residue-P4  phase  code  sequence  to  detect  the  target  at  the  maximum  unambiguous  range 
is  presented  (given  the  fixed  unambiguous  range  Ru  =21 50  m ).  The  relationship 

between  the  required  SNRRo  and  the  average  CW  power  Pcw  is  given  by  Equation  (26). 

In  Figure  65,  the  CW  signal  power  as  a  function  of  the  required  SNRRo  for 
mx  =  3,  m2  =  4,  and  m3  =5  is  plotted  with  BRj  =  1  MHz,  tb  =  1  /us,  t'h  =  0.333  ps, 
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G,  =  Gr  =  30  dB,  fc=  3  GHz,  A  =  0.1  m,  oy  =  100  m2,  FR=  5  dB,  and  L2  =  LRT  =  LRR  =  1, 
and  RnrdX  =  Ru  =  cMRSNSt'h  !  2  =  2150  m  (maximum  detection  range  without  ambiguity). 


Figure  65.  Average  power  of  the  CW  transmitter  for  robust  symmetrical  residue-P4 

with  ml  =  3,  m2  =  4,  and  m2  =  5. 

Note  that  the  difference  of  the  average  power  used  in  each  channel  is  due 
to  the  PG.  The  channel  m3  =  5,  giving  PG  =  2m  =  10,  requires  less  power  than  the  other 

channels,  m1  =  3  ( PG  =  6)  and  m2=  4  ( PG  =  8  ). 

c.  Unambiguous  Detection  Range  for  a  Constant  Output  SNR 

For  a  single  P4  code,  the  unambiguous  range  is  limited  to  only  one  code 
period  duration  Ru  =  cT  /  2  =  ctbNc  /  2.  In  this  section,  we  compare  the  maximum 

unambiguous  range  using  the  robust  symmetrical  residue-P4  phase  code  sequence  with 
the  corresponding  single  P4  code  sequences.  The  relationship  between  the  average  signal 
power  Pcw  and  the  maximum  detection  range  Rmy.  is  given  by  Equation  (27). 
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In  Figure  66,  the  detection  range  is  plotted  as  a  function  of  the  required 
CW  power  using  N  =  3  RSNS  for  ml  =  3,  m2  =  4,  and  m3  =  5  with  a  constant 

SNRRo=  13  dB,  B  =  1MHz,  tb=  l//s,  t'h  =  0.333  //s ,  G,  =  G,.  =  30  dB ,  /f=3GHz, 
A  =  0.1m,  <jt  =  100  m2 ,  T0  =  290  K ,  FR  =  5  dB ,  and  L2  =  LRT  =  =  1 .  Also,  the 

performance  of  each  P4  code  using  the  corresponding  Nc  is  plotted.  For  the  P4  code 
from  Equation  (28),  Ncl  =  NcRSNS  /  N  =  18/3  =  6,  Nc2  =  8,  and  Ae3=10.  The  single  P4 
codes  have  a  limited  unambiguous  range;  Ncl=  6  gives  Ru  =3x1 0s  x  1 0  6  x  6  /  2 
=  900  m,  Nc2  =  8  gives  Ru  =  3 x  10s  x  10~6  x8  /  2  =  1200  m,  and  Nci  =  10  gives  Ru  = 
3xl08  xlO"6  xlO/2  =  1500  m.  These  values  are  less  than  achieved  using  the  N  =  3 
robust  symmetrical  residue -P4  phase  code  sequences,  RU/i  =  cMRSNSt'b  /  2  =  2150  m. 


Figure  66.  Comparison  of  the  maximum  unambiguous  range  of  CW  radar  system  for 
SNRRo  =  13  dB  using  robust  symmetrical  residue-P4  code  for  m1  =  3, 
m2  =  4,  and  m3  =  5  with  the  radar  using  each  P4  code  individually  for  the 
corresponding  N  =  N  IN. 

1  G  c  CRSNS 
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In  conclusion,  the  robust  symmetrical-P4  phase  code  sequence  presented 
in  this  chapter  can  be  used  to  extend  the  unambiguous  range  of  CW  polyphase  radar 
systems.  The  issue  concerned  with  the  range  bin  error  due  to  the  fluctuation  of  the 
compression  outputs  is  also  resolved.  The  integer  Gray  code  properties  can  prevent  an 
error  in  the  detected  robust  symmetrical  residues.  In  the  next  chapter,  the  concluding 
results  from  using  the  number  theoretic  transforms  from  Chapters  III  to  V  to  extend  the 
unambiguous  range  in  CW  polyphase  radar  systems  are  discussed. 
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VI.  CONCLUDING  RESULTS  AND  RECOMMENDATIONS 


A.  EQUATION  SUMMARY 

The  summary  of  equations  used  to  compute  the  signal  and  range  detection 
parameters  for  each  modular  phase  code  (residue-Frank  phase  code,  symmetrical  residue- 
P4  phase  code,  and  robust  symmetrical  residue-P4  phase  code)  that  results  in  extending 
the  unambiguous  range  are  shown  in  the  Table  10.  The  equations  are  used  for  any  of  the 
coprime  system  moduli. 


Table  10.  Summary  of  equations  for  each  modular  phase  code. 


Parameters 

RNS-Frank 

SNS-P4 

RSNS-P4 

Phase  sequence 

(0k  ’0k') 

2n  k- 1 

0  i  = -  —  R  u 

Ym,k  m  1  m  ]  m>k 

-[-{S">*-2  JJ  4m 

ti,k'  =  ■f^RSm,k'-mf  ~\m 

Subcode 

period  (tb,t’b) 

h  ~  CPP !  fc 

h  =  CPP !  fc 

t'b=tbl  N  =  cpp/ Nfc 

Signal 

bandwidth  (B) 

Mtb 

1  ltb 

l/tb=l/(Nt'b) 

The  number  of 

subcodes 

(NC,NC  ) 

m2 

m 

2mN 

Processing  gain 

( PG ) 

m2 

m 

2  m 

Code  period  (T) 

2, 
m  tb 

mth 

2  niNt'h 
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Using  a  set  of  moduli,  we  can  extend  the  maximum  detection  ranges  beyond  those 
of  individual  moduli.  The  equations  for  calculating  the  dynamic  range,  the  number  of 
code  periods  within  the  dynamic  range,  the  range  resolution,  and  the  unambiguous  range 
are  shown  in  Table  1 1 . 


Table  1 1 .  Summary  of  equations  for  N  -  channel  modular  number  system. 


Parameters 

RNS-Frank 

SNS  -P4 

RSNS  -P4 

Dynamic  range 

(MrnS’  ^SNS’^RSNS^ 

N 

Wm 

i= 1 

L  J  N 

min]  n  mi,+  \  |  «/,  • 

[21=2  l=j+ 1  'J 

(even) 

3  2  15 

—  m,  H - m.  +  1  for 

2  1  2  1 

N  —  3  moduli  of  the 

[  1  J  1  N 

min  n  W//+  n  min 

l2 1=1  2l=j+^  J 

J-  (odd) 

form  2Z  -1,2Z,2Z  +1 

The  number  of 

M rns 

M RSNS 

code  periods 

m 

m 

2  Nm 

(Npj„) 

Range 

SL "l 

< 

resolution  (AR) 

2 

2 

2 

Unambiguous 

cMth 

cM-SNsh 

RSNsh 

range  ( Ru ) 

2 

2 

2 

In  the  SNS  and  RSNS,  the  dynamic  range  is  not  necessarily  greater  with  an 
increasing  value  of  the  moduli  values.  For  example,  the  dynamic  range  of  a  3-channel 

SNS  system  with  moduli  [7,  13,  16]  is  MSNS  =  69.  On  the  other  hand,  the  dynamic  range 
of  a  3-channel  RSNS  system  with  moduli  [7,  13,  19]  is  only  M RSNS  =55.  That  is,  ml  and 
m2  are  equal,  but  the  last  RSNS  modulus  is  /«,  =19,  which  is  larger  than  the  last  SNS 
modulus  m3  =16.  However,  the  dynamic  range  of  the  RSNS  is  less  than  the  SNS. 
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B.  COMPARISON  OF  3-CHANNEL  MODULAR  NUMBER  SYSTEMS  WITH 
THE  SAME  DYNAMIC  RANGE 

In  Table  12,  a  comparison  of  each  modular  number  system  having  the  same 
dynamic  range  (60),  the  same  subcode  period  (1  ps),  the  same  unambiguous  range  (9000 
m),  and  the  same  range  resolution  is  shown.  Note  the  bandwidth  for  the  RSNS-P4  system 
is  smaller  since  it  is  contains  N  phase  values  that  are  equal. 


Table  12.  Comparison  of  each  modular  phase  code  in  which  each  has  the  same 

dynamic  range  =  60,  and  tb,  t'b  =  1  jus  . 


Parameters 

RNS-Frank 

SNS-P4 

RSNS-P4 

Set:  Dynamic  range 

60 

60 

60 

Set:  Subcode  period  (ps) 

1 

1 

1 

Bandwidth  (MHz) 

1 

1 

0.33 

Unambiguous  range  (m) 

9,000 

9,000 

9,000 

Range  resolution  (m) 

150 

150 

150 

Moduli 

[3,  4,  5] 

[7,13,29] 

[3,  5,  7] 

Numbers  of  phase  values 

[3,  4,  5] 

[4,  7,15] 

[4,  6,  8] 

code  period  (ps) 

[9,16,  25] 

[7,13,29] 

[18,30,  42] 

Number  of  code  periods 

[6.67,3.75,2.4] 

[8.6,  4.6,  2.1] 

[3.3,2,  1.4] 

The  RNS  compared  to  the  SNS  is  more  efficient  in  covering  Ru  since  fewer  code 
periods  are  required  in  the  compression.  The  phase  values  from  the  phase  detector  are 
also  less  for  the  same  Ru .  Comparison  with  the  RSNS  shows  that  it  is  the  most  efficient  of 

all  schemes  due  to  the  fewer  number  of  phase  values  used,  the  smallest  number  of  code 
periods  required  in  the  compression  and  its  Gray  code  properties  that  limit  the  range  error 
to  one  range  bin.  In  terms  of  complexity,  the  size  of  the  compressor  is  longest  for  the 
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RSNS  since  the  dynamic  range  incorporates  the  fewest  code  period  (assuming  Nr  =1). 
The  SNS  requires  the  fewest  code  periods  in  the  compressor. 

C.  COMPARISON  OF  3-CHANNEL  MODULAR  NUMBER  SYSTEMS  WITH 
THE  SAME  MODULI 

In  Table  13,  each  number  theoretic  transform  that  has  the  same  set  of  moduli  [3  4 
5]  with  the  subcode  period  equal  to  1  ps  is  shown. 


Table  13.  Comparison  of  each  modular  number  system  in  which  each  has  the  same 

moduli  [3  4  5],  and  th,  t'b  =  1  jus 


Parameters 

RNS-Frank 

SNS-P4 

RSNS-P4 

Set:  Moduli 

[3  4  5] 

[3  4  5] 

[3  4  5] 

Dynamic  range  (units) 

60 

11 

43 

Set:  Subcode  period  (ps) 

1 

1 

1 

Bandwidth  (MHz) 

1 

1 

0.33 

Unambiguous  range  (m) 

9,000 

1,650 

6,450 

Range  resolution  (m) 

150 

150 

150 

Number  of  phase  values 

[3,  4,  5] 

[2,  3,  3] 

[4,  5,  6] 

Number  of  code  periods 

[6.67,  3.75,  2.4] 

[3.7,  2.8,  2.2] 

[2.4,  1.8,  1.4] 

In  the  RNS  phase  code,  the  extended  unambiguous  range  is  maximized  due  to  its 
maximum  dynamic  range.  In  the  SNS  phase  code,  the  number  of  phase  values  is  lowest 
due  to  the  construction  of  the  SNS  integers  within  the  modulus.  However,  it  gives  the 
shortest  unambiguous  range  since  its  dynamic  range  is  the  smallest.  The  RSNS  gives  the 
next  largest  unambiguous  range  of  the  three  schemes.  The  Gray  code  properties  make  it 
an  attractive  choice. 
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D.  RANGE  DETECTION  ERROR  CONTROL 


As  shown,  only  the  robust  symmetrical  residue -P4  phase  code  sequence  could 
prevent  a  range  error  due  to  the  noise  changing  the  phase  of  the  return  signal.  This  is  due 
to  the  Gray  code  property  that  the  RSNS  has.  For  the  other  modular  phase  codes  (residue  - 
Frank  and  symmetrical  residue-P4),  the  range  error  is  large  when  the  target  straddles  two 
range  bins.  The  higher  set  of  moduli  might  help  control  the  range  error  due  to  its  higher 
processing  gain.  However,  increasing  the  modular  number  would  increase  the  number  of 
phase  steps,  resulting  in  increasing  receiver’s  complexity  and  a  longer  code  period, 
resulting  in  higher  compression  loss  and  required  dwell  time  on  targets. 

E.  RECOMMENDATIONS  FOR  FUTURE  WORK 

The  modular  phase  codes  used  in  this  thesis  were  shown  to  be  able  to  extend  the 
unambiguous  range  in  CW  polyphase  modulation  radar  systems.  The  applications  of 
these  modular  phase  codes  in  different  types  of  radar  systems  would  be  interesting  to 
study.  Also,  the  performance  of  each  moduli  set  with  the  different  trade-offs  should  be 
further  considered.  The  Doppler  tolerant  property  that  the  P4  code  has  could  also  be 
investigated  for  the  symmetrical  residue-P4  and  the  robust  symmetrical  residue-P4  phase 
code  sequences.  Finally,  the  hardware  implementation  has  to  be  examined  and  a 
simulation  that  includes  hardware  imperfections  and  errors  must  be  conducted. 
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APPENDIX  -  LIST  OF  VARIABLES 


There  are  many  variables  used  throughout  this  thesis.  All  of  the  important 
variables  are  listed  here  to  provide  unambiguous  meanings. 


Table  14.  List  of  variables. 


Variable 

Meaning 

a 

Range  bin  index 

c 

Speed  of  light  (in) 

cpp 

Number  of  carrier  cycles  within  a  subcode 

fa 

Carrier  frequency 

fs 

Sampling  frequency  (Hz) 

h 

Integer  position  in  RSNS 

k 

Phase  index 

kB 

Boltzman’s  constant  (joule/K) 

m 

Coprime  modulus 

mi 

Coprime  moduli 

h 

Subcode  period  (s) 

B 

Bandwidth  (Hz) 

br, 

Receiver’s  input  bandwidth  (Hz) 

fr 

Receiver  noise  factor 

Gr 

Gain  of  receiving  antenna 

G, 

Gain  of  transmitting  antenna 

Grr 

Loss  between  the  radar’s  antenna  and  receiver 

Lrt 

Loss  between  the  radar’s  transmitter  and  antenna 

l2 

Two-way  atmospheric  transmission  loss 

M 

Frank  frequency  step 

^  RNS 

Dynamic  range  of  the  RNS 

Msns 

Dynamic  range  of  the  SNS 

^  RSNS 

Dynamic  range  of  the  RSNS 

N 

Number  of  coprime  modular  sequences 

K 

Number  of  subcodes  within  a  code  period 

N 

1 y cRSNS 

Number  of  subcodes  within  a  code  period  (RSNS) 

NP 

Number  of  code  periods  returned  from  target 
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Variable 

Meaning 

Np,m 

Number  of  code  periods  used  for  modulus  m 

K 

Number  of  receiver  reference  codes 

PG 

Processing  gain 

PSL 

Peak-to-sidelobe  level 

P 

1  RNS 

Period  of  a  single  RNS  sequence 

p 

1  SNS 

Period  of  a  single  SNS  sequence 

P 

1  RSNS 

Period  of  a  single  RSNS  sequence 

R 

Range  to  target  (m) 

K 

Unambiguous  range  (m) 

Ru,RNS 

Unambiguous  range  of  RNS  (m) 

R u, SNS 

Unambiguous  range  of  SNS  (m) 

R u.RSNS 

Unambiguous  range  of  RSNS  (m) 

^ m.k 

RNS  residue  kth  for  modulus  m 

RSm,k 

RSNS  symmetrical  residue  Mi  for  modulus  m 

c 

^  m,k 

SNS  symmetrical  residue  Mi  for  modulus  in 

SNR„ 

Input  signal-to-noise  ratio 

snrRo 

Output  signal-to-noise  ratio 

T 

Code  period  (s) 

T 

nij 

Code  period  for  coprime  modulus  mi 

V 

Target  velocity  (m/s) 

T 

Roundtrip  time  delay  (s) 

TpW 

Pulse  width 

A 

Signal  wavelength  (m) 

V 

Doppler  frequency  (Hz) 

A R 

Range  resolution  (m) 

(7  j 

Target’s  radar  cross  section  (m  ) 

0m,k 

Phase  code  for  modulus  m  with  index  k 
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